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SYMBOLS 
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uana 


A,  B,  C,  D ,  E  Coefficient  of  nth  order  polynomial  from  high  to  low.  As"  Bs"'^  ...Es^ 

b  Wing  span 

c  Wing  chord 

e  Mean  Aerodynamic  chord  (mac) 

Cd  D 

Drag  coefficient,  — 
qb 

Lift  coefficient, 

qs 

_  S£ 

Aerodynamic  non-dimensional,  rolling  moment  coefficient, 

C'e  Change  in  non-dimensional  rolling  moment  coefficient  with  change  in 

variable  i  (where  i  -  V,  p,  r,  &a.  &r).  Axis  system  is  other  than  principal  axis, 

C't  ■  ({Cft  +  Iw  C'„/I*x)/(1  -  Ixx'Ixx  Ia)l 

C  M 

Aerodynamic  non-dimensional,  pitching  moment  coefficient, 

C  N 

"  Aerodynamic  non-dimensional,  yawing  moment  coefficient, 

C',,  Change  in  non-dimensional  yawing  moment  coefficient  with  change  in 

variable  i  (where  i  -  V,  p,  r,  &a,  &n)-  Axis  system  is  other  than  prindpai  axis, 

C'ni  •  ((C,  +  Ija  C(i/lzz)/0~  Ixyixx  Ia)I 

P  Iz 

^''PovN  Dynamic  Departure  Susceptibility  Parameter,  Cn^  cosa  -  (r-)B  Ci^  sina 


Aerodynamic  non-dimensional,  side  force  coeffident,  — 

qb 
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'^x-'-af 

^  aCm 

IT 
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Cx.--^ 

p.  _  aCm 

Cx.--^ 

_  aCm 

*4v; 

Rolling  moment 

Yawing  momani 

SMaJotcfi 

Sf 

“-.uSb 

q,Sb 

ac, 

"■af* 

^  aCn 

""  \£L 
°2v; 

p  ac, 

'‘15: 

r  « 

""  I5I 

Cy-i2l 

^  aCn 

_  aCr 
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SYMBOL  DEFlNtTlQN  IMLS 

DIMENSIONAL  AERODYNAMIC  DERIVATIVES  DEFINED  BY  CHODY  (REFERENCE  (36)) 


yJ  =  (Yp  +  Do)  cospo  +  (Dp  -  Yo)sinpo  + 

( g  [cospoCcosoo  sinOo  -  sinoo  cosOq  oos^)  -  sinpo  cosBo  siiK^]] 

Vsec 

*  1 

Yp  =  sinoo  +  ^Ypcospo  +  DpsinPo) 

rad 

*  •! 

Yr  =  -  COSoto  +  ^YrCOSP  +  DrSinPo) 

rad 

Y*  =  ^(cospoCsinoo  cosOo  sin(t)o)  +  cospo  coseo  cos<|)o  ] 

Vsec 

Lp  =L'p 

Lp  «L'p 

Vsec 

Lp  =  L'p  +  (  +  QoJxy  )/^x 

Vsec 

Lf  =  L'r  +  (  -Poixy  +  Qo  ( ly  “  I2)  “  2RoIx2yix 

Vsec 

«  * 

Lfl  =  Lp  cosog  +  Lr  sinoo 

Vsec 

Np  =N'p 

Vsec* 

Np  =N'p 

Vsec 

N  J  =  N'r  +  (  Roly,  +  Qo  ( Ix  -  ly)  +  2PoIxyyiz 

Nr  =  N'r  +  (  -Qola  Fo^yz 

Vsec 

* 

Nq  -  N'p  cosoo  +  N'r  sinoo 

^-4ec 

<!>;  '1 

- 

(|>*  -  tanOo  cos<|)o 

rad 

(j)*  -  Po  tanOo 

rad 
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DIMENSIONAL  AERODYNAMIC  DERIVATIVES  DEFINED  BY  CHODY  (REFERENCE  (36)) 

Note,  if  L'p  and  N'ji  are  not  "small"  then, 

L'r  =  L'r  +  L'p  COSOo 
Up  =  Up  +  Up  sinoo 
N'p  =  N'p  +  N'p  sinoo 
N'r  =  N'r  +  N'p  COSOo 

C(„  Cm,  Dimensionless  body-axis  ( other  than  principal)  roll  and  yaw  moment  coeffi¬ 

cient 

D  Aerodynamic  Drag  lb 

F  Force  lb 

g  Gravitational  acceleration  constant,  32.174  •Vsec® 

h  Altitude  ft 

i  State  Variable,  i.e.,  a,  p,  6,  etc 

I  Identity  Matrix 

2 

Ix,  ly,  Iz  Body-axis  moments  of  inertia  slug-ft 

Ixz  Body-axis  Product  of  inertia  slug-ft^ 

jo)  Imaginary  part  of  Laplace  transform  variable  s  -  o  jco 

(yr  Distance  from  the  aircraft  c.g.  to  the  X-coordinate  of  the  vertical  tail  ft 

aerodynamic  center. 

L  Aerodynamic  Lift  lb 
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SYMBOL 
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Li.  Mi.  Ni 

(dmidi)  •  {1>W 

(a(MiVdi)  •  (H) 

(d(Niyai)(M.) 

1 

ft-  sec 

where  1  ■  v 

i-p.  a 

1 

sec^  -  rad 

i  =  p.  q.  r.  a.  p. 

1 

sec  -  rad 

i  =  6i 

1 

sec^  -  rad 

a'i 

Total  incremental  change  in  rolling  acceleration  due  to  incremental 
change  in  state  variable  quantity, 

(as  applicable) 

“  (Sfi  +  ^  Ni)/(1  - 

m 

Aircraft  mass 

slugs 

M 

Mach  number 

- 

M 

Pitching  moment  about  aircraft  body  axis 

ft-lb 

N 

Yawing  moment  about  aircraft  body  axis 

tt-lb 

N'i 

Total  incremental  change  in  yawing  acceleration  due  to  incremental 
change  in  state  variable  quantity. 

(as  applicable) 

-  (Ni  +  ^  u  y(i  - 

Ni 

Aircraft  transfer  function  numerator  relating  output,  i,  to  control  deflection,  j,  (i 

->  p,  q,  r,  etc;  j  &a.  &«.  M 

P 

Perturbed  roll  angular  velocity  about  the  X-body  axis 

racVsec 

P 

Total  roll  rate  about  the  X-body  axis,  P  -  Po  p 

rad/sec 

q 

Perturbed  pitch  angular  velocity  about  tire  y-body  axis 

racVsec 

Q 

Total  pitch  rate  about  the  y-body  axis,  Q  -  Qo  q 

radVsec 

q 

1  9 

Dynamic  pressure, 

lb 

ft2 
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SYMBOLS  (CONT) 

SYMBOL 

DEFINITION 

utuis 

r 

Perturbed,  yaw  angular  velocity  about  the  y-body  axis 

rad/sec 

R 

Total  yaw  rate  about  the  z-body  axis,  R  -  Ro  r 

racVsec 

s 

Laplace  operator  (s  ••  o  ijo}) 

" 

S 

Aircraft  wing  area 

ft2 

t 

Time 

sec 

T 

Thrust 

lb 

Tr 

First  order  time  constant  of  the  roll  subsidence  mode 

sec 

Ts 

First  order  time  constant  of  the  spiral  mode 

sec 

Te,,  Tea 

First  order  time,  constants  of  the  conventional  low  and  high  frequency  zeros  of 
the  pitch  attitude  numerator 

sec 

Tea 

First  order  time  constant  of  the  pitch  altitude  numerator  which  results  from  un- 
symmetrical  flight 

sec 

T^,,  T<()j 

First  order  time  constants  of  the  overdamped  roll  attitude  numerator 

sec 

U 

Control  vector 

- 

u 

Control  perturbation  vector 

- 

u 

Perturbed  linear  velocity  component  along  the  x-axis 

ft 

sec 

U 

Total  linear  velocity  along  the  (X  body  axis),  U  -  Uo  u 

ft 

sec 

V 

Perturbed  linear  velocity  component  along  the  (Y  body  axis) 

V 

Total  linear  velocity  along  the  (Y  body  axis),  V  -  Vo  v 

ft 

sec 

Vt 

Total  linear  velocity,  •/  +  V*  +  W® 

ft 

sec 

w 

Perturbed  linear  velocity  component  idong  the  (Z  body  axis) 

ft 

sec 

W 

Total  linear  velocity  along  the  (Z  body  axis),  W  Wo  w 

ft 

sec 

w 

Weight 

lb 

X 

State  vector 

. 
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SYMBOLS  (CONT) 

SYMBOL 

DEFINITION 

UNITS 

X 

State  perturbation  vector 

- 

X.  Y.Z 

Aerodynamic  forces  along  the  OXYZ  body-axis  system 

ft) 

Xi.  Yj.Zi 

I  d( )  /  di  ]/m  where  i  -  u,  v,  w. 

1 

where  i  -  p,  q,  r.  a,  P 

sec 

ft 

sec  -  rad 

where  i  &j 

ft 

sec^  -  rad 

GREEK 

DEFINITION 

UNUS 

a 

-1  W 

Angle  of  attack,  tan  — 

deg 

p 

-1  V 

Sideslip  angle,  sin  ^ 

deg 

y 

Flight  path  angle 

deg 

M 

Roll  angle  about  the  velocity  vector 

deg 

p 

Mass  density 

slugs 

ft® 

r 

Wing  geometric  dihedral  angle 

deg 

6i 

Surface  deflection  with  subscript 

deg 

A 

Transfer  function  denominator 

- 

Ci 

Damping  ratio  with  subscript 

- 

Eigenvalue 

- 

♦.e.v 

Conventional  perturbed  Euler  angles  (roll,  pitch,  yaw) 

deg 

e.  Y 

Conventional  Euler  angles  between  inertial  axis  and  aircraft  body  axis  (roll, 
pitch,  yaw) 

deg 

o 

Real  part  of  Laplace  transform  variable,  a  "•  a  t  jco 

racVsec 

B() 

Natural  frequency  of  the  denominator  or  numerator  root  with  subscript 

racVsec 

CO^ 

Rotational  rate  vector  of  reference  frame  A  with  respect  to  reference  frame  B 
coordinates 

radl/sec 

Sweepback  angle  at  the  quarter  chord 

deg 
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SYMBOLS  (CONT) 

MATHEMATICAL  DEFINITION 

_A_  Square  diagonal  Matrix  whose  diagonal  elements  are  distinct  eigenvalues 

A  define 

« 

^ ^  derivative  with  respect  to  time,  d/dt 

«>  infinity  (without  bound) 

d  Partial  derivative 

<  less  than 

>  greater  than 

II  II  norm 

( )  Denotes  Vector  quantity 

^  j  Matrix  Column  Vector 

( ]  Square  Matrix 

T  Matrix  T  ranspose 

Transfer  functions  are  presented  using  the  following  shorthand  notation, 
(s  +  a)  =  (a),  where  a  ■  ’A 

I^S^  +  ^ci)nS  +  (Onj  =  IC,  0)n] 
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SYMBOLS  (CONT) 

Subscripts  Description 


a  Aileron 


APR 

B 

COP 

DR 

DYN 

e 

I 

L 

( 

m 

n 

o 

P 

r 

R 

s 

SP 

SR 

VT 


Apparent 
Body  axes 

Coupled  stability  axis  system  (of  Kaiviste.) 
Dutch-roll 

Dynamic  stability  axis  system  of  Kaiviste 

Elevator 

Inertial  axes 

Left 

Aerodynamic  moment  about  the  x-axis 
Aerodynamic  moment  about  the  y-axis 
Aerodynamic  moment  about  the  z-axis 
Reference  or  nominal  value 
Phugoid 
Rudder 

Right,  Roll  mode 
Stability,  Spiral  mode 
Short  period 

Coupled  spiral-roll  subsidence  (lateral  phugoid) 
Vertical  Tail 
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SYMBOLS  (CONT) 

ACRONYMS 

DEFINITiON 

ARI 

Aileron-Rudder  Interconnect 

BVR 

Beyond  Visual  Range 

CIC 

Close-ln-Combat 

DOF 

Degree-of-Freedom 

OOC 

Out-of-Control 

STl 

Systems  Technology  Inc. 

WVR 

Within  Visual  Range 
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THE  BODY  AXIS  SYSTEM  IS  DEFINED  BY  THE  FOLLOWING: 


X0  longitudinal  body  axis  in  the  plane  of  symmetry  of  the  aircraft, 
positive  forward; 


Y  Q  lateral  body  axis  perpendicular  to  the  plane  of  symmetry  of 
the  aircraft,  usually  taken  in  the  plane  of  the  wings,  positive 
toward  the  right  wing  tip; 


Zg  vertical  body  axis  in  the  plane  of  symmetry  of  the  aircraft,  per¬ 
pendicular  to  the  longitudinal  and  lateral  axes,  positive  down. 


Figure  i(a)  Definition  of  Body  Axis  System  (Reference  (1)) 


XX 
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(6)  Stability  Axis 


Xs  longitudinal  atabili^/  axis,  paraim  to  the  projaction  of  tho  total 
voiodty  (Vt)  on  tha  plana  of  aymmatry  of  tha  aircraft,  positiva 
forward: 


Ys  latw’al  stability  axis,  coincidant  with  artd  positiva  in  tha  sama 
dlracUon  as  tha  latarai  body  axis; 


Z«  vortical  stability  axis  in  tha  plana  of  symmatry  of  tha  aircraft, 
^rpandicular  to  tha  longitudinal  and  latarai  axas,  positiva 
down. 


Figure  i(b)  Definition  of  Stability  Axis  System  (Rafaranca  (1)) 
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AERODYNAMIC  ANGLES 


a  pitch  angle  of  attack,  angle  b«tw«enth«Xraxis  and  ttwXB' 
axle,  positive  rotates  the  ■•>Z-axis  into  the  -fX-axis: 

a  >  tan‘^  (W/U),  -lao®  iai  180®; 


p  angle  of  sideslip,  angle  between  the  total  velocity  (Vt)  ertd  Its 
projectkxi  on  the  XZ-piane,  positiva  rotates  the  Wre^s  into 
the  ♦X.-axis:  p  -  sin*  ^  (VAAr).  -90®  i  p  S  90®; 


Figure  ii  Aerodynamic  Angles  (Reference(1 )) 
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Xb 

X| 


NOTE:  2i  IS  COINCIDENT  WITH  THE 
LOCAL  GRAVITY  VECTOR. 


naiPMTATlOW  ANGLES: 

y*w  angle,  positive  dockwise  about  the  4-z>axis  direction; 


pitch  angle,  positive  clockwise  about  the  ♦Y-aws  direction; 
roll  angle,  positive  dockwise  about  the  ♦X*axls  dkection; 

V,  0,  and  ♦  term  a  system  of  three  angles  which  uniquely  defines  Ihe  orleotationa  of  the  Body 

Axes  with  respect  to  an  Inertial  reference  set  d  axes.  Any  orientation  of  the  Body  Axis  System  Is  obtained 

(uniquely)  by  rotatlonally  displadng  It  from  the  referenoe  system  through  each  of  the  three  anjHes  in  turn. 

The  order  of  rotation  Is  Important  and  Is  defined  to  be  yaw-pitch-rol. 


Figure  ill  Earth-Body  Orientation  (Euler)  Angles  (Reference(1 )) 
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INTRODUCTION 


BACKGROUND 

Out'of-Control^  (OOC)  flight  is  not  new  to  the  Navy  fleet.  Aircraft  departure^  (from  controlled  flight)  and  OOC 
accidents  have  plagued  aviation  from  if  s  infancy.  The  importance  of  being  able  to  confidently  predict  when 
an  aircraft  has  the  potential  to  depart  controlled  flight  is  vividly  illustrated  using  the  F- 14  as  an  example.  From 
1973  through  1987,  thirty  F*14  aircraft  have  been  lost  due  to  OOC  accidents.  This  equates  to  one*  third  of 
all  F-1 4  accidents.  Not  only  is  aircraft  departure  resistance  an  important  safety  issue,  but  it  is  also  paramount 
in  maintaining  mission  effectiveness  in  today’s  close-in-combat  (CIC)  environment.  The  CIC  environment 
today,  and  in  the  near  future,  is  projected  to  be  characterized  by:  (1)  short  duration  maneuvering,  i.e.,  time 
compression,  (2)  an  increase  in  importance  placed  on  the  transition  between  beyond-visual-range  (BVR) 
and  within-visual-range  (WVR)  maneuvering  and  (3)  expansion  of  the  combat  arena  into  the  low-speed 
post-stall  flight  envelope  regime.  These  three  characteristics  require  harder  maneuvering  and  improved 
agility  for  today’s  and  future  aircraft.  Improved  agility  over  today's  current  fighters  will  resiJt  in  aircraft  capable 
of  higher  turn  rates,  increased  acceleration  and  deceleration  capability  and  probably,  most  importantly,  skew 
maneuvering  capability  involving  nose  pointing  at  the  expense  of  energy  conservation.  To  make  these 
improvements  in  maneuvering  possible  requires  that  the  aircraft  possess  some  degree  of  departure 
resistance.  Understanding,  predicting  and  designing  for  departure  resistant  aircraft  (without  sacrificing 
desired  levels  of  maneuverability)  has  spurred  the  development  of  many  aircraft  d^>arture  susceptibility 
criteria  and  guidelines  that  originate  with  the  development  of  the  basic  directional  weathercock  stability 
derivative,  Cn«-  This  report  will  expand  on  each  of  the  major  departure  susceptibility  criteria,  and  highfight 
each  of  their  applications  and  imitations.  A  more  detailed  aircraft  departure  susceptible  analysis  approach 
is  recommended  to  assess  the  effects  of  asymmetric  flight  and  maneuvering  dynamics. 


1  Sm  Appendix  A,  “Qlossary  Of  Defined  Terms* 
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OBJECTIVE 

The  primary  objective  of  the  study  reported  herein  is  to  provide  guidelines  for  the  analytical  determination  of 
aircraft  departure  susceptibility  applicable  to  the  designs  of  future  fighter  aircraft. 

SCOPE 

This  report  entails  a  brief  review  of  the  derivation,  application  and  implied  limitations  of  the  aircraft  departure 
susceptibility  criteria  presented  in  the  MIL-STD-1 797A  Flying  Qualities  specification  (Reference  (2)).  To  begin 
the  discussion,  the  rigid  body  six  degree-of-freedom  equations-of-motion  and  aircraft  stability  concepts  are 
reviewed,  focusing  on  their  application  to  determine  the  susceptibility  of  an  aircraft  to  depart  from  controlled 
flight.  Finally,  a  methodology  is  suggested  that  allows  a  designer  to  determine  the  departure  susceptibility  of 
an  aircraft  configuration  through  the  course  of  its  design  (up  to  and  including  flight  test)  with  increasing  levels 
of  accuracy  as  the  design  progresses. 

METHODOLOGY 

In  providing  guidelines  for  the  analytical  determination  of  aircraft  departure  susceptibility,  the  intent  was  not 
to  endorse  one  criteria  over  another,  or  to  rehash  what  MIL-STD-1 797A  recommends  or  even  to  propose  a 
"new  and  improved"  criteria.  Rather  an  attempt  is  made  to  point  out  the  departure  susceptibility  prediction 
techniques  which  are  appropriate  at  various  stages  of  the  design  phase  as  a  function  of  the  required  accuracy. 
As  pointed  out  in  Reference  (3),  consistent  simple  prediction  techniques  are  important  early  in  the  design 
phase  to  establish  the  configuration  and  make  trade  studies  between  performance  and  departure  resistance. 
Later  in  the  development  cycle,  increased  accuracy  is  required  to  avoid  costly  configuration  changes  and 
reduce  risk  during  flight  test. 

There  are  four  primary  contributors  of  aircraft  departure  from  controlled  flighL  They  are:  (1)  high  angle-of-at- 
tack  bare  airframe  aerodynamic  flight  characteristics,  (2)  use  of  aerodynamic  flight  controls  (i.e.,  pilot  in  the 
loop),  (3)  inertial  coupling  and  (4)  kinematic  coupling  (see  Figure  1).  The  influence  of  each  of  these  factors 
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I  Aircraft  Bara  Airframe  Hlgh>AOA  AarodyiMunle  Flight  Charactariatics 

•  Nonlinear  «vith  raapact  to  angla-of-attaok  and  tidaalip 

•  Flow  breakdown  and  Adverra  Vortex  Shedding  Effects  Common 

Associated  Causes  of  Departure 

1.  Aircraft  unstable  rSrecdonally  with  stable  dihedral  effect  (or  vice  versa); 
Aircraft  may  depart  but  is  not  likely  to  be  divergent 

2.  Aircraft  is  unstable  directionidly  and  has  unstable  dihedral  effect 
divergent  departure  likely 

II  Use  of  Aerodynamic  Flight  Contras 

•  Aileron  and  rudder  effectiveness  greatly  reduced 

Associated  Causes  of  Departure 

1 .  Use  of  aileron  may  aggravate  situation  due  to  adverse  yaw  generated  at  high 
angle^*attack  becoming  the  dominant  control  effect 

2.  Use  of  prolonged  or  misapplication  (cross>oontrolling)  of  control  inputs  at  high 
angle-^-attack  could  induce  departure. 

III  Inertial  Coupling 

•  Moments  generated  due  to  inertial  coupling  becoming  more  pronounced  at  large 

angles-of-attack. 

Associated  Cause  of  Departure 

1 .  Inaease  in  inertial  coupling  effect  at  large  angle-of-attack  Ikely  to  place  the 
aircraft  in  a  flight  condition  more  susoepUbie  to  departure. 

IV  Kinematic  Coupling 

•  Sideslip  angle  (angler-attack)  generated  due  to  kinematic  coupling  can  become 

more  pronounced  with  inaeasing  angle-of-attack  (sideslip). 

Associated  Cause  of  Departure 

1 .  From  a  departure  susoeptit^lity  viewpoint  any  generation  of  large  amount  of 
sideslip  (or  angle-of-attack)  is  undesirable  because  it  has  the  potential  to  place 
the  aircraft  in  a  flight  condition  that  is  more  susceptible  to  departure. 

Figure  1  Primary  Causes  of  Airaaft  Departure  from  Controlied  Right  (Reference  (4)) 
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as  they  contribute  to  the  susceptibility  of  an  aircraft  to  depart  controlled  flight  has  been  investigated  by  many 
researchers  to  date.  In  addition,  Kalviste  has  recently  addressed  each  of  the  miyor  departure  causing 
elements  with  the  exception  of  closed-loop  pilot/flight  control  effects  in  a  single  methodology  (See  reference 
(5)). 

RESULTS  AND  DISCUSSION 

General 

To  achieve  the  accuracy  at  which  to  predict/analyze  aircraft  stability  parameters  and  flying  qualities  at  high 
angles-of-attack  depends  chiefly  on  the  ability  to  formulate  a  'Valid*  aerodynamic  mathematical  model.  This 
then  assumes  that  the  desired  quantities  (i.e.,  Ci^,  C^,  etc)  can  either  be  accurately  measured  in  a  wind  tunnel 
or  be  calculated  analytically  with  a  good  confidence  level. 


Of  particular  interest  to  the  flight  dynamicist's  in  terms  of  predicting  an  aircraft’s  tendency  to  depart  from 
controlled  flight  are  those  flow  phenomena  that  significantly  vary  with  angle-of-attack  and  that  can  cause 
asymmetric  effects  even  when  the  aircraft  maintains  a  zero  sideslip  attitude.  Two  of  the  most  important 
phenomena  of  this  kind  are  (1)  the  formation  and  asymmetric  shedding  of  forebody  vortices  and  (2)  the 
formation  and  asymmetric  bursting  of  wing  leading-edge  vortices.  These  phenomena  become  even  more 
complex  if  an  oscillatory  motion  is  superimposed  on  the  primary  steady  flight  trajectory  (Reference  (6)). 
These  high  angle-of-attack  flow  phenomena  have  large  effects  on  all  of  the  aerodynamic  characteristics 
of  the  aircraft  including,  of  course,  the  static  and  dynamic  stability  parameters.  The  most  important 
of  these  effects  on  the  dynamic  stability  parameters  are  (1)  large  nonlinear  variations  of  stability 
parameters  with  angle-of-attack,  angle-of-sideslip  and  rate  of  coning  as  well  as  with  amplitude  and 
frequency  of  oscillation,  (2)  significant  aerodynamic  cross-coupling  between  longitudinal  and  lateral 
degrees  of  freedom  (i.e.,  Cm;,  CnJ,  (3)  time-dependent  and  hysteresis  effects  and  (4)  strong 
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configuration  dependence.  Reference  (6)  discusses  each  of  these  four  effects  in  great  length,  but  for  the 
purposes  of  this  report,  the  focus  will  be  on  the  impact  each  of  these  four  elements  have  on  the  methodology 
required  to  accurately  predict  aircraft  departure  susceptibility. 

To  illustrate  both  the  magnitude  and  the  suddenness  of  the  nonlinear  variations  in  dynamic  derivatives  with 
angle-of-attack,  the  yaw  and  roll  damping  derivatives  for  the  X-31 A  aircraft  are  presented  in  Figure  2.  The 
yaw  damping  derivative  (Figure  2a)  exhibits  a  very  sudden  and  very  large  (twice  the  low-a  value  and  the  sign 
reversed)  unstable  peak  at  an  angle-of-  attack  of  40  degrees  while  the  roll  damping  derivative  (Figure  2b) 
exhibits  an  equally  large  and  sudden  variation  (again  about  twice  the  low-a  value  and  the  sign  reversed)  with 
angle-of-attack  such  that  it  is  undamped  over  the  28  to  50  degree  AOA  range.  Work  done  at  NASA  Langley 
(References  (8)  and  (9))  suggest  that  the  primary  mechanism  for  these  effects  is  associated  with  the  flow 
phenomena  emanating  from  the  aircraft's  forebody.  The  studies  showed  the  peak  instabilities  to  be  largely 
independent  of  the  wing  sweep  angle  and  the  presence  of  vertical  tails. 

In  addition  to  dynamic  derivatives  exhibiting  large  nonlinear  variations  with  angle-of-attack,  they  have  also 
been  shown  to  exhibit  significant  nonlinear  variations  with  the  frequency  and  amplitude  of  oscillation.  This 
variation  takes  place,  for  the  most  part,  where  the  derivatives  change  suddenly  with  angle-of-attack.  An 
example  of  the  variation  of  the  roll  damping  derivative  with  oscillation  amplitude  is  shown  in  Figure  3  for  a 
fighter  configuration  at  low  subsonic  speed  (Reference  (10)).  The  large  unstable  peak  that  occurs  at  an 
angle-of-attack  of  35  degrees  when  the  amplitude  of  oscillation  is  1 5  degrees  decreases  at  larger  amplitudes 
and  completely  disappears  at  an  amplitude  of  ±  20  degrees.  With  the  knowledge  of  these  frequency  and 
amplitude  effects,  NASA  Langley  presented  the  forced  oscillation  test  data  of  the  X-31  A  aircraft  with  an 
envelope  that  encompassed  the  variability  of  the  data  due  to  the  amplitude  and  frequency  variations  during 
the  forced  oscillation  testing.  Note  that  these  are  not  "bands  of  uncertainty"  but  envelopes  due  to  actual  data 
taken. 

Figure  4  illustrates  the  range  of  data  for  the  yaw  and  roll  damping  derivatives  for  the  X-31  A  aircraft. 
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Figure  2a  Yaw  Damping  Derivative  Versus  Angle-of-Attack  for  the  X>31  A  Aircraft  (Reference  (7)) 
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Figure  2b  Roll  Damping  Derivative  Versus  Angle-of-AttacK  for  the  X>31  A  Aircraft  (Reference  (7)) 
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Figure  3  Effect  of  Oscillation  Amplitude  on  Damping  in  Roll  (Reference  (10)) 
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Figure  4  Range  of  Data  for  the  Yaw  and  Roll  Damping  Derivatives  of  the  X-31 A  Aircraft  (Reference  (7)) 
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Additionally,  the  nominal  and  range  of  data  of  the  cross-coupling  dynamic  derivatives,  Cnp  and  C^,  as  a 
function  of  angle-of-attack,  are  shown  in  Figures  5  and  6  respectively. 

In  cases  where  the  amplitude  and  frequency  of  oscillation  effects  are  significant,  the  derivative  concept 
becomes,  at  best,  a  "good"  guess  of  the  dependence  of  aerodynamic  damping  with  angle-of-attack. 

Aerodynamic  Cross  Coupling 

The  introduction  of  cross-coupling  derivatives  requires  simultaneous  consideration  of  the  complete  six 
degrees-of-freedom  equations  of  motion  of  an  aircraft,  rather  than  the  decoupled  sets  of  equations  for  the 
longitudinal  and  lateral/directional  degrees-of-freedom.  The  phenomena  of  aerodynamic  cross-coupling  can 
occur  at  high  angles-of-attack  when,  1 )  lateral  aerodynamic  reactions  (such  as  caused  by  the  lateral  motions 
of  the  forebody  vortices)  may  occur  on  an  aircraft  as  a  result  of  some  longitudinal  motion  such  as  pitching 
or  vertical  translation  or  2)  longitudinal  aarodynamic  raactions  (such  as  those  caused  by  the  longitudinai 
motion  of  the  vortex  burst  locations)  may  occur  on  an  aurcraft  as  a  result  of  some  lateral  motion  such  as 
rolling,  yawing  or  lateral  translation.  Static  cross-coupling  data  for  an  early  version  of  the  X-31 A  aurcraft  is 
presented  in  Figure  7.  Figure  7a  illustrates  the  pitching  moment  as  a  function  of  amgle-of-attack  and  sideslip 
angle,  while  Figure  7b  illustrates  the  yawing  moment  coefficient  as  a  function  of  angle-of-attack  and  sideslip 
angle.  For  this  version  of  the  X-31  A  aircraft,  it  is  apparent  that  aircraft  aerodynamic  pitching  moment  is  a 
strong  function  of  sideslip  angle  above  -15  degrees  (i.e.,  Cmp)-  Likewise  the  aircraft  aerodynamic  yawing 
moment  is  a  strong  function  of  angle-of-attack  (i  «-.  CnJ.  This  aerodynamic  cross-coupling  phenomena  is 
important  and  must  be  modelled  to  correctly  analyze  aircraft  stability  characteristics  at  high  angles-of-attack. 
For  further  discussion  on  aerodynamic  cross-coupling  see  reference  (1 2).  In  reference  (1 3)  an  analysis  was 
conducted  to  determine  the  influence  of  the  Mp  term  on  the  aircraft  dynamics  of  the  F-4  aircraft.  This  was 
done  by  calculating  the  e/6e,  (>/5a  and  r/&r  transfer  functions  from  the  six  degree-of-freedom  coupled 
(P  *  0)  equations-of-motion.  For  this  particular  case  (F-4  Aircraft;  oo  -  20°;  Po  "  1° ).  Mp  principally  affects 
the  lateral  and  longitudinal  short-period  modes  with  a  resultant  increase  in  short  period  damping  (Csp)  and 
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YAW  DUE  TO  ROLL  RATE,  1/rad 


ANGLE  OF  ATTACK,  Alpha  -  deg 


Figure  5  Yawing  Moment  Coefficient  Due  to  Roll  Rate  for  the  X-31 A  Aircraft  (Nominal  and  Range  of  Data) 
(Reference  (7)) 
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Figure  6  Rolling  Moment  Coefficient  Due  to  Yaw  Roll  Rate  for  the  X-31 A  Aircraft  (Nominal  and  Range  of 
Data)  (Reference  (7)) 
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Figure  7  (Concluded)  Static  Aerodynamic  Cross-Coupling  Data  for  the  Preliminary  Configuration  of  the  X-31  Aircraft  (I 
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a  decrease  in  Dutch-roll  damping  (Cdr)  (sea  Rgure  8).  Further  it  appears  that  with  Mp  -  -1 .63,  closure  of  any 
or  all  of  the  three  loops  shown  in  Figure  8  could  easily  result  in  an  unstable  Dutch-roll  mode,  (Reference 
(1 3)).  For  Mp  positive,  the  migration  of  the  longitudinal  short-period  and  Dutch-roll  modes  due  to  the  infhiertoe 
of  Mp  is  reversed  from  that  shown  in  figure  8.  That  is  the  longitudinal  short-period  moves  toward  the  right 
(unstable)  while  the  Dutch-roll  moves  to  the  left. 

The  concept  of  aerodynamic  cross-coupling  (Ci,,  Cm,)  with  respect  to  dynamic  derivatives  has  only  recently 
been  introduced  (See  Reference  (14)).  But  on  the  basis  of  sensitivity  studies  as  described  in  references  (1 5) 
and  (16),  it  has  been  shown  that  the  inclusion  in  the  equations  of  motion  of  the  dynamic  cross-coupling 
derivative  C^,  can  cause  instabilities  in  the  Dutch-roll  and  roll  subsidence  modes  of  motion.  As  with  any 
second  order  effect,  it  must  be  kept  in  mind  that  the  effect  of  the  cross-coupling  derivatives  is  quite  dependent 
on  the  remaining  stability  characteristics  of  the  aircraft.  The  smaller  the  static  margin  and  the  lower  the 
aerodynamic  damping,  sudi  as  represented  by  Cip,  C^,  Cn,,  Cn^  derivatives,  the  ntore  sensitive  the  aircraft 
motion  will  be  to  variations  in  cross-coupling  derivatives  and  vice  versa. 

Time-Dependent  Effects 

In  addition  to  modelling  quasi-steady  effects  such  as  represented  by  dynamic  derivatives  due  to  angular 
velocities  (Cip,  Cn,.  Cpp  etc.),  the  existence  and  importance  of  modelling  the  purely  unsteady  effects  (as 
represented  by  derivatives  due  to  the  time  rate  of  change  of  angular  deflections,  C^,  Cn^,  etc)  at  high 
angles-of-attack  is  documented  in  references  (17)  and  (18).  The  translational  derivatives  as  they  are  often 
referred  to,  constitute  part  of  the  dynamic  results  obtained  with  standard  wind-tunnel  techniques  of  oscillating 
a  model  around  a  fixed-axis.  These  tests  always  gives  composite  derivatives  expressed  as  (Cn,-  Cn^  cosa) 
and  the  like  (see  Table  I).  Common  practice  when  formulating  this  data  into  the  standard  Taylor  series 

expansion  of  the  aerodynamic  forces  and  moments  is  to  ignore  the  a  and  p  effects  (or  to  introduce  a  simple 
correction  for  them)  and  to  use  the  composite  derivatives  in  place  of  the  purely  rotary  ones.  At  low 
angles-of-attack  the  errors  introduced  by  this  formulation  is  often  small  and  the  simplification  is  large  enough 

to  make  it  justifiable.  At  higher  angles  of  attack,  however,  the  a  and  p  effects  can  become  quite  substantial 
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Table  I  Dynamic  Moment  Derivatives  (Reference  (6)) 
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and  may  no  longer  be  ignored  or  corrected  for  in  a  simple  fashion  (Reference  (6)).  A  more  detailed  discussion 
on  the  modelling  and  significance  of  the  time-dependent  effects  can  be  found  in  reference  (19). 

Based  on  a  series  of  experiments  conducted  in  the  Tracor  Hydronautics  Ship  Model  Basin,  the  pure 
longitudinal  and  lateral  rotational  and  translational  derivatives  were  measured  for  a  0.19  scale  model  of  a 
preliminary  X-31 A  configuration  operating  at  a  Reynolds  number  of  about  1 .28  x  10^  From  Reference  (20) 
the  following  results  are  cited. 

1 .  The  measured  value  of  the  composite  pitch  damping  derivative  (Cnv,  Cma)  Is 

equal  to  the  algebraic  summation  of  the  individual  rotational  (Cnvq)  and  trans¬ 
lational  (CmJ  up  to  24  degrees  angle-of-attack.  Beyond  24  degrees  angle-of- 
attack  there  is  a  marked  difference  between  the  results.  At  70-degrees 
angle-of-attack  the  algebraic  sum  of  Cnv,  >s  about  twice  the  value  of 

the  measured  composite  derivative  derived  from  the  combined  pitching  mo¬ 
tion.  This  is  shown  in  figure  9. 

2.  Figure  1 0  and  1 1  show  the  measured  values  of  the  composite  derivative 
(Cvr  -  Cyp)s  and  (Cc,  -  Ctp)s  compare  quite  well  with  the  algebraic  sum  of  the 
individual  coefficients  for  the  entire  angle-of-attack  range. 

3.  In  the  case  of  the  composite  yaw  damping  derivative  (Cn, -Cnp)s  (see  Figure 
12),  agreement  with  the  algebraic  sum  of  (Civ)5  and  (Cn^s  is  good  up  to  an 
angle-of-attack  of  approximately  30  degrees.  For  angles-of-attack  greater 
than  30  degrees  the  derivative  derived  from  the  algebraic  sum  of  (Ctv)s  and 
(Cnp)s  has  a  much  smaller  positive  (unstable)  value  than  the  values  obtained 
from  the  combined  yawing  motion.  It  is  apparent  that,  depending  on  how  the 
aircraft  aerodynamics  is  mathematically  modelled,  the  predicted  behavior  of 
the  aircraft  dynamic  response  in  the  lateral  mode  could  be  substantially  dif¬ 
ferent. 
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ANCLE  OF  ATTACK  a  IN  DECREES 


Figure  9  Comparison  of  the  Composite  Pitch  Damping  Derivative  (Cnv,  +  Cma )  and  the  Sum  of  the  Pure 
Pitching  (Cm, )  and  Pure  Plunging  (Cmi )  Derivatives  with  Angle-of-Attack  for  a  Preliminary 
Configuration  of  the  X-31 A  Aircraft  (Reference  (20) 
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ANCLE  OF  ATTACK  a  IN  DECREES 


Figure  10  Comparison  of  the  Composite  Sideforce  Damping  Derivative  (Cy,  ~  and  the  Sum  of  the 
Pure  Yawing  (Cy,)8  and  Pure  Sideslipping  (Cy^H  Derivatives  with  Angle-of-Attack  for  a  Prelimi¬ 
nary  Configuration  of  the  X-31 A  Aircraft  (Reference  (20)) 
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ANCLE  OF  ATTACK  a  IN  DECREES 


Figure  1 1  Comparison  of  the  Composite  Derivative  (C^  -  C^)s  and  the  Sum  of  the  Pure  Yawing  (C^)s  and 
Pure  Sideslipping  (C(^)s  Derivatives  with  Angle-of-Attack  for  a  Preliminary  Configuration  of  the 
X-31 A  Aircraft  (Reference  (20)) 
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Figure  12  Comparison  of  the  Composite  Yaw  Dcunping  Derivative  (Cn,  -  Cn0)s  and  the  Sum  of  the  Pure 
Yawing  (Cn,)*  and  Pure  Sideslipping  (Cn«)t  Derivatives  with  Angle-of-Attack  for  a  Preliminary 
Version  of  the  X-31 A  Aircraft  (Reference  (20)) 
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Hysteresis  Effects 

High  angle-of-attack  flow  phenomena  such  as  asymmetric  vortex  shedding,  vortex  burst,  or  periodic 
separation  and  reattachment  of  the  flow  are  frequently  responsible  for  aerodynamic  hysteresis  effects.  Such 
hysteresis  is  characterized  by  a  double*valued  behavior  of  the  steady-state  aerodynamic  response  to 
variations  in  one  of  the  motion  variables  such  as  angie-of-attack,  angle  of  sideslip  or  spin  rate.  In  the  presence 
of  hysteresis,  the  dynamic  derivatives  measured  in  large  amplitude  oscillation  experiments  may  have  two 
distinct  components.  Namely,  one  associated  with  the  small-amplitude  oscillation,  and  a  second  one 
representing  the  effect  of  the  hysteresis. 

Configuration  Dependence 

The  intricate  vortex  pattern  that  exists  around  an  aircraft  configuration  at  high  artgles-of-attack  is  very  sensitive 
to  even  small  changes  in  aircraft  geometry,  (see  References  (3),  (21 )  and  (22)).  The  forebody  vortices  are  greaHy 
dependent  on  the  planform  and  the  cross-sectional  geometry  of  the  aircraft  nose,  as  well  as  on  the  presence  of 
any  protuberances  on  the  forebody  (such  as  nose  boom,  strakes,  etc.)  that  may  affect  the  stability  of  an  existing 
vortex  pattern.  Such  protuberances  give  rise  to  new  vortices  and  create  and/or  alter  oorKlitions  of  vortex 
interactions. 

The  wing  leading-edge  vortices,  in  addition  to  being  a  strong  function  of  the  leading-edge  sweep,  are  also 
known  to  be  greatly  affected  by  various  modifications  of  the  wing  itself,  such  as  the  addition  of  wing  fences. 
All  of  these  variations  of  the  geometry  of  the  wing  affect  not  only  the  position  and  strength  of  tire  wing  vortices, 
but  also  the  all-important  location  at  which  these  vortices  breakdown.  Knowledge  of  the  sensitivity  of  the 
forebody  and  wing  vortex  systems  to  even  small  geometric  changes,  and  in  turn  the  strong  dominance  of 
aircraft  lateral-directional  stability  at  high  angles  of  attack  on  forebody  aerodynamics,  is  indicative  of  how 
carefully  one  must  approach  analyzing  aircraft  open-loop  stability  at  high  angles  of  attack.  In  reference  (6), 
a  preliminary  assessment  was  made  of  the  relative  significance  of  modelling  each  of  the  dynamic  derivatives 
for  an  aircraft  at  high  angles-of-attack.  The  results  of  the  assessment  are  given  in  Table  II. 
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Table  II  Relative  Significance  of  Dynamic  Moment  Derivatives  at  High  Angles-of-Attack  (Reference  (6)) 


Type  of  derivatlve(a) 

Derivativefs) 

Aircraft 

Significant 

Miaailes 

Direct 

Cmq.  Cof,  Cip 

Yas 

Yas 

Cross 

Cnpi  Cl, 

Yas 

Yes 

Cross'ooupling 

C|q,  Coq 

Yes 

Yes 

Cross'ooupling 

Cmp 

No 

Yes 

Cross*ooupling 

Cnif 

No 

No 

Aooeleration 

Cm^i  Ci^,  C|^ 

Yas 

? 

Acceleration 

Cn^ 

? 

? 

NADC-00048^ 


STABILITY  AND  THE  RIGID  AIRCRAFT  EQUATIONS-OF-MOTION 

Before  discussing  analytical  methods/criteria  for  predicting  aircraft  departure  susceptibility,  an  up  front 
definition  of  departure  from  controlled  flight  is  essential.  The  term  "departure"  is  generally  defined  to  be  a 
divergent,  large  amplitude,  uncommanded  aircraft  motion  (i-e..  pitch>up,  nose-slice,  roll  reversal,  etc.).  This 
definition  can  be  further  refined  and  a  distinction  be  made  between  open-loop  and  closed-loop  departures. 
Open-loop  departures  refer  to  departures  from  controlled  flight  that  are  due  to  instabilities  in  the  basic  aircraft. 
That  is,  even  if  the  pilot  does  not  move  the  controls,  small  perturbations  in  the  aircraft  states  build-up  until 
the  aircraft  can  no  longer  be  controlled.  In  a  closed-loop  departure,  the  basic  aircraft  may  or  may  not  be 
unstable,  but  the  addition  of  a  pilot-loop  closure  creates  an  unstable  vehide/pilot  system.  With  the  exception 
of  the  lateral  control  departure  parameter  (LOOP)  and  l/T^i,  the  departure  susceptibility  criteria  presented 
in  MIL-STD 1797A,  solely  addressed  open-loop  instabilities  of  the  basic  airframe  dynamics. 

Analytically  predicting  the  occurrence  of  an  aircraft  departure  from  controlled  flight  is  invariably  based  on  the 
assumptions  of  the  equations-of-motion  (EOM)/math  model  used  to  describe  the  response  of  the  aircraft  to 
initial  conditions,  pilot  control  inputs,  atmospheric  disturbances  and  system  failures  (i.e.,  propulsion,  flight 
control  system,  etc.).  As  can  already  be  inferred  by  9ie  definition  presented  for  an  open-loop  departure, 
departure  from  controlled  flight  connotes  an  instability  with  respect  to  an  initial  flight  condition.  While  a 
closed-loop  departure  has  to  do  with  the  stability  of  the  aircraft  due  to  pilot  control  inputs  and  loop  closures. 
Like  all  physical  systems,  the  six  degree-of-freedom  EOM  describing  the  motion  of  a  rigid  aircraft  are 
nonlinear  and  have  time  varying  parameters  to  some  degree  (see  Figure  13).  Associated  with  nonlinear 
systems,  the  concept  of  stability  is  very  much  different  than  that  of  linear  systems.  For  a  linear  system, 
stability,  or  the  lack  of  it,  is  a  property  of  the  linear  constant  parameter  system,  and  the  stability  is  determined 
by  the  eigenvalues  (solution  of  the  homogenous  EOM)  of  the  system.  In  addition,  the  stability  of  a  linear 
system  is  unaffected  by  the  initiai  conditions  or  the  forcing  functions  (care  must  be  taken  when  applying  this  last 
statement  that  one  does  not  violate  assumptions  that  permitted  the  use  of  a  linear  model  in  the  first  place). 
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U-^^RV-QW-gsine 

(1) 

V  ■  ^  ♦  PW  -  RU  ♦  goose  sin^ 

(2) 

W  -  ^  ^  QU  -  PV  >  goose  oos^ 

(3) 

W 

Q-^lM-(p2-R2}l«-RP(I*-Ix)l 

(5) 

(6) 

V  -  (Q  aln^  *  R  oos^)seoe 

(7) 

e  -  Q  oos^  -  R  sin^ 

(8) 

^  -  P  ■»  Q  sin^tane  ♦  R  oos^tarte 

(9) 

Note:  The  force  (X.Y,Z)  and  moment  (S(,M,N)  terms  represent  the  external  forces  and  moments  acting  on 
the  vehicle  other  than  gravity  (l-e.,  aerodynamic,  engine,  etc.) 

Figure  13  Six  Degree-of-Freedom,  Nonlinear,  Rigid  Body  Equations-of-Motion  (Body-Axis;  lyz  ■  Ixz  ■  0) 
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In  contrast,  a  nonlinear  system  may  be  stable  for  one  input  or  initial  condition  and  unstable  for  another.  That 
is,  stability  for  a  nonlinear  system  is  generally  not  a  system  property  but  a  property  of  a  particular  solution 
corresponding  to  a  particular  initial  condition  of  the  nonlinear  system.  For  this  reason,  there  unfortunately 
exists  no  convenient  method  for  the  determination  of  the  stability  of  a  nonlinear  system  in  general  (the  Routh, 
Horwitz,  and  Nyquist  stability  criteria  are  only  applicable  to  linear  systems).  In  most  cases  this  forces  attention 
to  focus  on  the  stability  of  a  nonlinear  system  in  the  'neighborhood*  of  a  solution.  The  reformulation  of  the 
nonlinear  system  equations  into  a  'linearized*  set  of  equations  about  a  defined  reference  point  is  an  effective 
and  widely  used  technique  to  begin  to  understand  the  stability  characteristics  of  a  nonlinear  system.  This 
technique  is  justifiable  provided  that  the  magnitudes  of  the  perturbation  quantities  and  the  chosen  operating 
point  used  to  define  the  linear  model  reproduces  the  nonlinear  system  dynamics  'accurately*  over  a 
"suHicient"  time  span.  When  the  linearizing  assumptions  are  not  Justified  the  analysis  of  the  mathematical 
model  must  be  carried  out  by  one  or  more  of  the  nonlinear  analysis  methods  (see  Reference  (23)).  An 
exception  to  this  exists  if  the  analyst  has  the  capability  of  deriving  multiple  linear  math  models  from  a  nonlinear 
system  math  model  using  a  computerized  numerical  method  technique.  Each  of  the  models  would  then 
provide  stability  information  valid  for  a  particular  region  of  the  state  space  of  interest. 

A  technique,  known  as  the  linear  approximation  theroms  of  Lyapunov,^  states  that  for  suffidentiy  small 
disturbances,  the  stability  (or  instability)  of  an  equilibrium  solution  of  a  nonlinear  system  can  be  determined 
in  all  cases  from  the  linear  model  if  none  of  the  eigenvalves  is  zero.  If  one  of  the  eigenvalves  is  zero  than 
the  set  of  equations  is  said  to  be  critical  according  to  Lyapunov  and  the  stability  of  the  nonlinear  solution  is 
dictated  by  the  discarded  nonlinear  terms. 

Note  that  this  concept  of  stability  addresses  the  initial  states  and  not  the  control  inputs.  Thus  this  definition  is 
applicable  to  investigating  the  potential  of  an  open-loop  departure  occurring.  On  the  other  hand,  the  stability 


'The  reader  may  wish  to  review  the  formal  definitions  and  theroms  of  the  small  disturbance  stability  theory 
provided  in  References  (24)  and  (25). 
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concept  in  terms  of  the  control  inputs  (and  not  the  initial  conditions)  is  described  by  the  foliowing  oortoept;  A 
linear  system  (described  by(:(t-Ax-»Bu:y-Cx)  whose  Laplaoe  transform  of  the  output  y(s)  is, 
y(s)  •>  C(SI-A)‘^  X  (o)  C(Sl-A)*^Bu(s) 


has  bounded-input  bounded-output  (BIBO)  stability  if 


II  u(t)  II  < «» 
X(0)-0 


lly(t)li<~ 


It  turns  out  that  time-invariant  linear  systems  are  BIBO  stable  if  all  eigenvalves  of  [A]  lie  in  the  open  left  half 
plane.  (Note  if  the  system  modes  are  all  completely  observable  and  controllable  then  Lyapunov  and  BIBO 
stability  are  one  in  the  same.)  In  some  cases  the  actual  response  of  y(t)  may  not  be  BIBO  stable  H  for  example 
an  excited  observable  mode  does  not  appear  in  the  transfer  function  (C(SI-A)'^B)  due  to  pole-zero 
canoellation  implying  the  unoontrollabiiity  of  a  normal  mode. 


The  existence  of  a  valid  'linearized”  model  of  a  nonlinear  system  however  is  not  always  assured.  Linearized 
equations  have  meaning  only  if  the  higher  order  terms  (products  and  cross  product  terms  of  the  perturbation 
variables)  are  "smaller”  than  the  linear  terms,  otherwise  the  higher  order  terms  can  dominant  the  system 
behavior. 


As  explained  earlier,  the  application  of  the  concept  of  stabiiity  (in  terms  of  Lypanov)  to  a  nonlinear  system 
described  by  equations  (1)  to  (9)  requires  studying  the  behavior  of  the  nonlinear  system  in  the  vicinity  of  a 
particular  sdudon.  The  particular  solution  of  interest  is  the  equilibrium  solution(s)  of  the  nordinear  system  for 
reasons  still  to  be  discussed.  Consider  now  reformulating  the  nonlinear  equation,  given  by  equations  (1)  to 
(9)  using  small-pertubation  theory  (See  figure  14).  To  do  this,  the  state  and  input  variables  are  redefined  in 
terms  of  the  sum  of  a  reference  state  (and  input)  quantity  (subscripted  ”0”)  arKi  a  perturbed  state  (and  input) 
quantity  Oower  case).  The  result  is  given  by  equation  (13)  of  figure  14,  where  A  and  B  are  assumed  to  be 
time-invariant  matrices  and  the  column  matrix  h(x)  contains  all  the  nonlinear  terms  as  a  function  of  the 
perturbed  states,  x. 
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ivwtth* 


ygM) 


R«d«fin«  th«  state  and  Input  variftbiM  such  that. 


X-  xo  +  x 


u  -  uo  u 


Equation  (10)  baoomas, 


x-/{x,u.t) 


y-9(x.u.t) 


or  in  steta-spaca  form  aquation  (12)  can  ba  wrtttan, 


X  -  Ax  Bu  4-  h(x) 


y-Cx 


whera:  A.  B  ara  ascumad  to  ba  tima-invafiant  matrioas  dafinad  by 


XDk  UO  Ma  I 


h(x)  -  noninaarterma  as  a  function  of  ttw  parturbad  states,  x. 


Figure  14  Ganaral  Formulation  of  the  Linear  State  Equations  for  a  Nonlinear  System 


NADC-90048-60 


According  to  the  small  disturbance  stability  theory  of  Lypanov,  the  stability  of  the  solution  x(t)  can  now  be 
determined  if  the  eigenvalues  of  [A]  are  nonzero  (positive  or  negative).  This  statement  holds  true 
provided  the  two  conditions  below  are  satisfied.  (Reference  (25)) 
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The  second  condition  above  implies  that  the  higher  order  nonlinear  terms  of  h(x)  are  negligible  compared  to 
the  linear  terms  of  A.  If  an  eigenvalue  of  the  derived  linear  system  is  zero  then  the  stability  characteristics 
must  be  determined  from  the  nonlinear  equations. 

There  are  two  reasons  why  the  reference  or  nominal  state  (xo)  used  in  small  perturbation  theory  is  chosen 

to  be  a  trimmed  condition  in  which  all  the  acceleration  terms  are  zero  (i.e.,  V-  0,  *  0).  The  first  is  that  the 

linear  system  equations  most  accurately  represent  the  nonlinear  system  dynamics.  This  is  shown  mathe¬ 
matically  by  equations  (14)  through  (17). 


Vb-(Vb)o 


+  AV 
A  ^ 


Now  if  ~  ”  0  is  satisfied, 
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The  second  reason  for  choosing  the  reference  state  as  a  trim  state  is  that  the  total  state  and  control  trajectories 
equate  to  the  sum  of  the  constant  nominal  values  and  the  linear  perturbation  time  history  values. 


^B(t)  =  (^B)o  +  A7B(t) 

(18a) 

waft)  =  (IJejo  +  A^Wt) 

(18b) 

U(t)  :  (U)q  +  AU(t) 

At  this  point,  the  nonlinear  rigid  aircraft  EOM  of  Figure  (13)  have  been  reformulated  in  the  perturbed 
state-space  form  of  equation  (1 3)  (see  Figure  1 5).  The  matrix  [A]  multiplied  into  [x]  along  with  the  [B]  matrix 
multiplied  into  [u]  represents  the  linear  portion  of  the  equations  and  the  column  matrix  h(x)  represents  the 
nonlinear  portion.  The  (A]  matrix  is  commonly  referred  to  as  the  stability  or  Jacobian  matrix  while  the  [B] 
matrix  is  referred  to  as  the  control  matrix. 

For  sufficiently  small  disturbances  the  product  and  cross  products  of  the  perturbation  variables  may  be 
considered  negligible  and  the  nonlinear  portion  of  the  equation,  h(x)  is  eliminated  leaving  only  the  matrix 
equation  that  describes  the  six  degree-of-freedom  linear  dynamics  for  a  general  reference  condition  (recall, 

that  this  method  is  utilized  based  on  the  premise  lim  »  0  (Reference  (25)). 

x-rO  II  ^  II 

The  (A]  matrix  given  in  figure  15  is  completely  general  with  the  only  assumptions  being  that  the  aircraft  is 
symmetrical  about  its  xz-plane  and  that  attitude  (6)  and  btmk  angle  (^)  perturbations  are  approximately  less 
than  1 5  degrees. 

Simplifying  assumptions  associated  with  the  complete  6  DOF  linear  EOM  are  used  in  many  flight  dynamic 
applications  to  reduce  the  complexity  of  the  problem  and  provide  engineering  insight.  The  problem  of 
predicting  aircraft  departure  from  controlled  flight  is  no  exception  as  will  be  seen  in  the  development  of  many 
of  the  existing  departure  susceptibility  criteria  parameters.  The  major  point  of  the  discussion  that  follows  is  to 
point  out  when  these  simplifying  assumptions  may  not  be  justified  and,  in  other  cases,  where  they  are  valid. 
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Rgure  15  Nonlinear  Rigid  Aircraft  Equations  of  Motion  in  Perturbed  State-Space  Form  (Body-Axis) 
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Figure  15  (Concluded)  Nonlinear  Rigid  Aircraft  Equations  of  Motion  In  Perturtod  State-Space  Form  (Body-Axis) 
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CHOICE  OF  THE  LINEARIZED  REFERENCE  FLIGHT  CONDITION 

Although  the  A-matrix  of  Figure  1 5  is  linear,  it  is  still  highly  complex  due  to  the  generality  presumed  for  the 
reference  trim  condition.  Because  of  their  complexity,  application  of  the  equations  in  their  complete  form  is 
seldom  used.  Instead,  simpler  trim  cases  that  still  reveal  stability  and  control  problems  are  more  commonly 
used.  Such  is  the  case  in  the  derivation  of  the  departure  susceptibility  parameters  CnpovN'  LCDP  .  a-p,  and 

5a.  The  derivation  of  these  parameters  assumes  steady  (^,io  -0),  straight  (y, >0),  symmetric 
( Vo.  Vo  0),  wings  level  «  0)  flight.  There  are  not  many  (analysts,  pilots,  etc.)  who  would  consider  steady, 
straight,  symmetric,  wings  level  flight  to  be  important  initial  conditions  practically  associated  with  departure 
susceptible  flight.  However  the  fact  remains  that  CpovN  LCDP  are  well  accepted  departure  susceptibility 
parameters  that  have  correlated  well  with  wind  tunnel,  and  flight  test  data  over  the  last  two  decades  (see 
References  (26)  to  (28)). 

There  are  significant  simplificaticns  and  insight  to  be  gained  if  a  reference  flight  condition  (Xo)  is  chosen 
such  that  Vo  (Po),  Po,  Ro  and  (|>o  are  zero.  Namely,  the  longitudinal  (u,  w  (or  a),  q  and  6)  and  lateral/directional 
(v(or  p),  p,  r,  linear  EOM  sets  can  be  decoupled  from  the  complete  linear  six  degree-of-freedom  aircraft 
EOM  (as  given  by  the  [A]  matrix  in  Figure  15)  and  solved  independently  of  each  other  (see  Figures  16  and 
17).  However,  as  shown  by  the  authors  of  References  (24)  and  (29)  through  (33),  these  simplifications  in 
many  cases  may  net  be  iustif  ied  and  often  their  solution  leads  to  misleading  stability  information.  In  the  worse 
case,  an  assumed  stable  flight  condition  may  in  fact  be  unstable. 

To  highlight  how  the  major  individual  elements  of  the  stability  matrix  are  influenced  by  asymmetric  flight 
conditions  (P,  0)  and  nonzero  angular  rates  (Po,  Qo,  Ro  *  0),  the  Stability  Matrix  ([A])  given  in  Figure  15 

has  been  simplified  and  presented  in  Figure  18.  Three  simplifying  assumptions  were  made.  They  are,  (1) 

the  cross  product  of  inertia,  (1xz)b,  is  zero,  (2)  the  vtP)ond\^a)  dynamic  translational  derivatives 
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Jacobian  (or  Stability )  matrix  lef/axl  evaluated  at  x 
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Figure  17  Decoupled  Lateral/Directional  Body-Axis  Stability  Matrix 
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Figure  18  Simplified  Linear  Stability  Matrix  Applicable  to  Asymetric  Maneuvering  Right 
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(S£v.  Nv .  Sfw>  Mw,  Nw,  Yv,  Zw)  are  zero,  and  (3)  the  aerodynamic  static  and  dynamic  longitudinal  •  lateraVdirec- 
tional  cross-coupling  derivatives  (S£w.  Mv,  Nw,  Nq)  are  zero. 

A  very  thorough  discussion  on  the  effects  of  altitude,  h,  velocity,  Vt.  the  aerodynamic  angles,  oo,  po 
asymmetric  flight  ^po  Oj  and  steady  maneuvering  (Pq,  Qo.  Po  0)  dynamics  on  the  unaugmented  rigid 
aircraft  stability  in  terms  of  the  eigenvalues  and  the  modal  shapes  is  given  by  Stengel  in  References  (30) 
and  (32).  The  general  findings^  of  this  work  are  summarized  in  Table  III. 


^  Note,  these  results  were  derived  for  a  reference  aircraft  (small  supersonic  fighter  type  designed  for  air 
superiority  missions)  but  can  be  generalized  in  most  cases. 
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Table  m  Summary  Of  The  Effects  Of  Altitude,  Velocity,  Aerodynamic  Angles  And  Steady  Maneuvering  On  The  Unaugmented  Rigid 
Aircraft  Stability.  (Reference  (30)  and  (32)) 
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Tr  is  changed  only  slightly 


Table  m  (Coat)  Summaiy  of  the  Effects  of  Altitude.  Velocity.  Aerodynamic  Angles  and  Steady  Maneuvering  on  the  Unaugmented  Rigid 
Aircraft  Stability.  (References  (30)  and  (32)) 
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through  u  and  w 


Table  m  (Cent)  Summary  of  the  Effects  of  Altitude.  Velocity.  Aerodynamic  Angles  and  Steady  Maneuvering  on  the  Unaugmented  Rigid 
Aircraft  Stability.  (References  (30)  and  (32)) 
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Table  m  (Concluded)  Summary  of  the  Effects  of  Altitude,  Velocity,  Aerodynamic  Angles  and  Steady  Maneuvering  on  the  Unaugmented 
Rigid  Aircraft  Stability.  (References  (30)  and  (32)) 
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CHRONOLOGICAL  REVIEW  OF  THE  DEVELOPMENTS  OF  AIRCRAFT  DEPARTURE  SUSCEP¬ 
TIBILITY  CRITERIA 

Much  has  been  written  on  the  subject  of  criteria  which  predict  aircraft  departure  susceptibility.  A  review  of 
these  criteria  immediately  points  out  that  linearization  of  the  nonlinear  rigid  body  aircraft  equations-of-motion 
(see  Figure  13)  is  common  to  the  derivation  of  each,  with  the  exception  of  Bihrie’s  Design  Plots. 

Table  IV  chronologically  cites  the  major  departure  susceptibility  criteria  developments  along  with  the 
associated  researcher(s)  and  the  referenced  work. 

It  should  be  noted  that  some  of  the  departure  susceptibility  criteria  developments  shown  in  Table  IV  are  more 
closely  classified  as  "stability"  parameters  (i.e.,  Cnacop>  opposed  to  parameters  that  have 

been  correlated  with  aircraft  departure  susceptibility. 

This  section  presents  a  review  of  the  departure  susceptability  criteria  developements  given  in  Table  IV.  Each 
of  the  criteria  will  be  discussed  ( in  chronological  order )  in  terms  of  its  basic  derivation,  the  impact  the  criterion 
has  on  aircraft  design  and  the  advantages  and  shortcomings  of  the  criterion. 


43 


Table  IV.  Chronology  of  the  Major  Aircraft  Departure  Susceptibility  Criteria  Developements 


NADC-90048-60 


DESCRIPTION/ 

COMMENTS 

Lateral/directional  body-axis  stability 
derivative  rsquirements  for  positive 
static  stabMty. 

Open-loop,  lateral-diiectional  require¬ 
ments  for  positive  static  stabNly  about 
the  velocity  vector. 

Closed-loop,  latoral-ditectlonal  require¬ 
ments  for  positive  static  stabNty  with 
respect  to  rol  oontroL 

ARemate  form  of  the  Cn,o^and  LCDP 
departure  parameters. 

•  o-S  >  0  implies  Cn„^  >  0 

•  a-  p  >  a»  implies  ICDP  >  0 

Criterion  derived  based  on  the  summa¬ 
tion  of  yawing  moments  about  the  z- 
stabilHy  axis  being  positive. 

First  criterion  to  oonelale  types  and 
aavarih#  of  aircraft  daoaitura  susoat}- 

tbil^  With  a  cross-plot  Of  the  open-loop 
latoraVdirectional  departure  parameter 
Cn^^  and  the  closed-loop  parameter 
LCDP. 

RESEARCHER(S)/ 

REFERENCES 

Unknown 

References:  24, 34 

Moul  &  Paulson  (NASA, 
Langley)  Referetioes:  3, 25, 
30,35, 36,37,38,39,40.41 

Jenny  (McDonnel 

Douglas)  References:  37. 
42.43 

WebSTTum  (Wlright- 

Paitarann  Air  Pnmt  Baaal 

Refeteinces:  27.  X,  37, 44, 
45,46 

Titiriga/Show  (Northrop) 
References:  47 

Johnston  (STI)  References: 
48 
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Table  IV  (Cent).  Chronology  of  the  Major  Aircraft  Departure  Susceptibility  Criteria  Developements 
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Table  IV  (Concluded).  Chronology  of  the  Major  Airaafl  Departure  Susceptibility  Criteria  Developements 
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(UNKNOWN):  Cnp,Cc,, 

PATE/RESEARCHER : 

<  1950;  (UNKNOWN):  Cnp.  Cjp 


CRITERIA : 

Positive  lateral/directional  aerodynamic  static  stability. 
Cnp  >  0 :  Basic  static  directional  weathercock  stability 
C(p  <  0 :  Dihedral  Effect  (lateral  static  stability) 


BASIS  OF  CRITERIA: 

Stable  static  lateral/directional  moment  data 


TYPE  OF  CRITERIA: 

Open-loop  Static  Stability  Criteria 


CRITERION  PLOT: 


DESIGN  IMPACT: 


Cnp :  •  Primarily  determines  natural  frequency  of  the  Dutch-roll  mode 

-  factor  in  determining  spiral  stability 

-  prevents  excessive  sideslip  and/or  yawing  motion  while  maneuvering  and/or 
in  turbulent  air. 

-  aids  in  turn  coordination 

Cfp :  -Prime  factor  in  damping  of  the  Dutch-roll  and  spiral  mode. 

•  Good  Cor  C(p  °  small  negative  (stable)  values 
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•  Stable  Spiral  Stability  ^  •  large  negative  values 

•  involved  in  maneuvering  characteristics  of  the  airframe  especially  with 
rudder  cJone  ojntrols  near  stall. 


The  sketches  below  illustrate  the  effect  of  varying  Cnp  and  on  the  lateral/directional  charac¬ 
teristic  modes. 


\  DUTCH 

Vmode 
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Figure  19  Lateral/Directional  Root  Locus  Plot  Varying  C(p  (Reference  (33)) 


IMAGINARY 
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DUTCH  ROLL  \ 
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NEGATIVE 

SPIRAL  MODE\ 

4. - —  real  axis 


Figure  20  Lateral/Directional  Root  Locus  Plot  Varying  Cn^  (Reference  (34)) 
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CRITERION  ADVANTAGE(S> : 

Data  is  available  early  to  affect  preliminary  design  phase  (i-®-*  “  /(^R'  -^**«  •  wing  loca¬ 

tion,  Cl,  aeroelasticity,  wing-fuselage  interference;  Cn^  -  /(Syr.  fuselage/wing  contribution). 

CRITERION  SHORTCOMINGS : 

Criterion  address  only  aerodynamic  static  stability  of  the  unaugmented  aircraft  with  respect  to 
sideslip  disturbances.  The  influence  of  maneuvering  dynamics  and  aerodynamic  cross-coupling 
are  not  considered. 
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MOUL  &  PAULSON:  CnpovN 
DATE/BESEARCHER  : 

1958;  Moul  &  Paulson  (NASA,  Langley) 


CRITERIA : 

A.  Positive  lateral/directional  open-loop  aircraft  static  stability  about  the  velocity  vector 

where : 


^'*i>ovN  *  CnpCOSa  “  |lxjg  (^®) 

B.  Aileron  Alone  Departure  Parameter  (AADP)/Lateral  Control  Departure  Parameter  (LCDP) 
Positive  closed-ioop  lateral-directional  stability  when  lateral  control  is  used  to  provide  constant 
roll  rate. 


where : 


LCDP(orAADP)>0 


LCDP- 


(Note,  LCDP  is  defined  with  stability-axis  derivatives.) 


(20) 


BASIS  OF  ^"PDYN  CRITERIA : 

Derived  from  the  general  aircraft  rigid  body  equations-of-motion  by  considering  the  stability  of 
the  linear  uncoupled  lateral/directional  dynamics  (i.e.,^  0  degrees  and  Vo  ~  0  ff/sec)  and 
assuming  the  reference  equilibrium  state  has  zero  angular  rates  (i.e.,  Po  *  Qo  ■■  Ro  ~  0  deg/sec 

and  yo  ~  ^  ”40  ~  0  deg/sec). 

To  test  for  stability,  Routh's  stability  criteria  is  applied  to  the  characteristic  equation,  ~  A)l  ”  0 
(where  [A]  is  given  in  figure  17,  See  equation  (21a)). 

AX^+BX®  +  CA.2+DA,+  E-0  (21a) 
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A  dynamic  system  becomes  unstable  and  divergence  will  occur  when  one  or  more  of  the  roots 
of  the  characteristic  equation  (eigenvalues)  becomes  positive.  According  to  Routh's  Stability 
Criterion,  the  characteristic  equation,  described  by  the  quartic  above,  will  have  positive  eigen¬ 
values  if  any  of  the  coefficients  (A,  B,  C.  D,  E)  or  if  the  combination  of  coefficients  BCD-AD^-B^E 
(Routh's  Discriminant)  becomes  negative. 

Solving  for  the  coefficients  of  the  characteristic  equation  yields, 

A-1 

B  ™  -(Yv  L'p  +  N'r) 

C  -  N'v  (Uo-Y,).L'v  (Yp  +  Wo)  +  L'p  (Yv  +  N'r)  +  (YvN'r  N'pL'r) 

D-NVlU(Yp  +  Wo)-YvL'pJ 

+  N'vlL'p(Yr-Uo)-L'f(Yp  +  Wo)l 
+  N'p  [YvL'r  -  L'v  (Yr  -  Uo)) 

+  g  COS0  [N'v  tanOo  -  L'v] 

E  -  -gcose  {-L'vN'p  tanOo  +  N'v  [Lp  tanSo  -  Lr )-  L'vN'r )  (21  e) 


(21b) 
(21c) 
(21  d) 


Moul  &  Paulson  concluded  that  the  oscillatory  instability  that  they  observed  was  associated  with 
a  change  in  the  sign  of  the  "C  coefficient  of  the  characteristic  equation. 
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Rearranging  and  simplifying  the  ”C”  coefficient  above  to  nondimensional  form  yields, 

C  -  CV«»«  ■  cvsira  ♦  ^  •  Cy,)  (22) 

+  (C't  C'n,  -  C',  C'opl 

The  "C  coefficient  given  in  equation  (22)  is  rewritten  as  shown  in  equation  (23)  to  aid  in  a  later 
discussion  on  recent  proposed  extensions  of  the  CnpQy^  concept  as  well  as  to  help  understand 
its  limitations. 

C  -  A  +  B  +  C  (23) 

where  :  A  =  C'npfcosa  -  ^  Cy,)-  •  C'tp  (sina  +  •  Cy,) 

C-CYplC'n,-^j{^j^C'^,l 

To  arrive  at  the  expression  for  CnpoY^  most  often  seen  in  the  literature  (e.g.,  Mil  Standard- 
1797A)  the  terms  of  equation  (23)  containing  products  of  derivatives  are  assumed  to  be  small 
compared  to  the  other  terms  (i.e.  this  drops  terms  B  and  C  from  equation  (23))  and  the  terms 

^^CYr  and  ^~CYp  are  assumed  small  in  ttie  A  term  compared  to  cosa  and  sina  respec- 

y  V  . 

tively. 
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Dropping  these  terms  from  equation  (23)  yields  the  expression, 


C-CneCOSa-N^ 


(g)^C^,sinc 


(24) 


Equation  (24)  is  referred  to  as  Chhiyn  because  it  has  been  found  to  be  a  good  indicator  of  ab¬ 
solute  (not  relative)  lateral/directional  dynamic  stability  (without  directly  calculating  the  aircraft 
eigenvalves,  i.e.,  Cnp^y^  <  0  implies  lateral/directional  instability)  for  a  host  of  aircraft  (see  Refer¬ 
ence  (26)). 


TYPE  OF  CRITERIA: 

Open-loop  Dynamic  Stability  Criteria 


CRITERION-PLOT : 


(Voea) 

1+) 

(-» 

DESIGN  IMPACT: 

The  CnpQy^  stability  parameter  shows  that  the  effective  directional  stiffness  is  a  function  not  only 
of  the  weathercock  stability  Cnp,  but  also  a  function  of  the  dihedral,  C^,  and  the  product  of  iner¬ 
tias  Ix,  Iz.  and  Ixz.  The  parameter  illustrates  clearly  how  at  high  angle-of-attack,  a  negative  Cnp 
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value  (possibly  due  to  tail  immersion  in  the  wing/body  wake)  can  be  stabilized  by  overriding 
negative  (stable)  values  of  (Iz/Ix  Ci,). 

As  an  approximation  only  (it  is  not  very  good  for  most  swept  wing  fighter  aircraft  which  have  rela¬ 
tively  large  values  of  and  significant  Ixz  values)  the  CnH^rN  ^  related  to  the  Dutch- 

roll  natural  frequency  as  shown  in  equations  (25)  to  (27)  below. 

anon  =  N'p,  +  YvN'r.  (25) 

(26) 

Provided  YvN'r, «  N'p„  then 

(2^ 


CRITERION  ADVANTAGE(S) ; 

Data  is  available  early  to  affect  preliminary  design  phase  (i.e.,  C|^  °;  /  ( AR^  V4>  >  wing  loca¬ 
tion,  etc ) ;  Cnp  =  /(SvT.  e^T.  «fc.)) 

CRIT£RLQN.SH0R.TCQM1NGS : 

The  ability  of  CnpovN  fo  predict  directional  divergence  can  be  affected  in  several  ways,  all  of 
which  are  important  to  understanding  its  application  and  limitations.  First  off,  it  must  be 
reiterated  that  the  CnpovN  parameter  is  derived  purely  from  application  of  the  Routh  Stability 
Criteria  to  the  "C"  coefficient  of  the  linear  decoupled  lateral/directional  equations-of-motion. 
Other  types  of  instabilities  can  also  occur  and  have  been  correlated  with  departure  suscep¬ 
tibility.  For  example,  equation  (27)  reveals  that  Cnpo^M  directly  relates  to  the  static  stability  of  the 

Dutch-roll  mode.  However,  Cnpoy^  >  0  does  not  predude  the  existence  of  dynamic  instability 
due  to  the  Dutch  roll  damping  term  being  negative.  This  could  be  revealed  by  analyzing  the 
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"B-coefficient"  of  the  characteristic  equation  (see  equation  (21a)).  Still  other  instabilities  as¬ 
sociated  with  longitudinal  and  lateral/directional  dynamic  coupling  have  been  associated  with 
departure  susceptibility.  The  properties  of  these  instabilities  will  be  discussed  further  under  the 
Chody  criteria  and  the  1  /Tea  departure  susceptibility  criteria. 

The  factors  that  directly  affect  the  accuracy  of  the  correlation  between  the  Cnpcxm  parameter  and 
departure  susceptibility  (i.e.,  lateral/directional  instability/divergence)  are: 

(1)  The  terms  that  were  assumed  negligible  in  reducing  the  C-coefficient  of  the 
linear  lateral/directional  equations-of-motion  to 

In  reference  (25),  calculations  were  made  for  two  different  configurations  to  com¬ 
pare  the  values  of  the  C-coefficient  (equation  (22))  and  the  Cnpoy^  parameter 
(equation  (24))  to  support  the  simplification.  In  terms  of  supporting  the 
simplifications,  the  reported  results  were  excellent  for  the  two  configurations 
considered.  That  is,  the  contribution  of  the  B  and  C  terms  of  equation  (23) 
were  negligible  compared  to  the  terms  of  A  (see  figure  21).  Additionally  the 
Cvp  and  Cvr  terms  of  the  A  term  are  also  negligible.  It  is  interesting  to  note 
that  the  comparisons  of  the  C-  coefficient  (unsimplified)  and  the  simplified 
CnpQY^  parameter  reproduced  in  figure  21  from  reference  (26)  was  con¬ 
sidered  up  to  an  angle-  of-attack  of  only  35-degrees.  In  comparison  today's 
high  performance  fighter  aircraft  are  capable  of  trimming  up  to  angles-of-at- 
tack  of  55  degrees  (e.g.,  McDonnell  Douglas  F/A-18)  and  experimental 
aircraft  such  as  the  X-31 A  are  being  designed  to  maneuver  up  to  angles-of-at- 
tack  as  large  as  70  degrees  to  achieve  maximum  performance  goals.  The 
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a.  Configuration  A 


Figure  21  Comparison  of  the  Lateraf/Directional  Characteristic  Quartic  C-Coefficient  with  Cn^gy^  for  Two 
Confjguratons  (Reference  (26)). 
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same  comparison  of  the  C-coefficient  with  the  CnpQY,^  parameter  has  been 
repeated  for  a  preliminary  design  of  the  X-31 A  aircraft  from  20  degrees 
up  to  an  angle-of-attack  of  70  degrees.  As  illustrated  in  Figure  22,  the 
X-31  A  data  comparing  the  values  of  the  C-coefficient  and  the  Cd^iyn 
parameter  overlay  each  other.  In  the  case  of  the  X-31  A  configuration 
the  static  characteristics  (  Ci^ ,  Cnp )  entirely  washout  the  second  order 
effects  of  the  dynamic  terms.  This  data  further  supports  use  of  the 
simplifying  equation  described  by  the  Cd^^yn  parameter  to  approximate 
the  C-coefficient  (equation  (22))  of  the  iateral/directional  quartic. 


As  pointed  out  in  Reference  (36),  the  fact  that  the  dynamic  derivatives  can  be¬ 
come  important  in  determining  the  C-coefficient  values  of  certain  aircraft  con¬ 
figurations  might  explain  the  inconsistency  of  the  empirically  derived 
boundaries  for  a  minimum  CnpoY^  value.  As  an  example,  the  work  of  Refer¬ 
ence  (47)  suggested  that  the  design  guidelines  of  Figure  23  be  used  in  lieu 
of  CnpoYN  ^  0  w^an  considering  the  presence  of  possible  aircraft  asym¬ 
metries,  destabilizing  external  loads  and  nonlinear  inertial  coupling  moments 
during  maneuvering  flight  (Reference  (36)).  Reference  (38)  simplifies  this 
even  further  and  suggests  that  Cnp^YN  ^  0.004(1 /deg). 


(2)  The  assumption  that  Ci  and  Cn  are  linear  with  sideslip  angle. 
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Key.  □  ;  C-Coeffident  of  the  Linear  Lateral/ 

Directional  Characteriatic  Quartic  (equation  (22^). 

O'  ^"♦oYN 


(1/DEG) 


0.01 


30  40  60  60 

Angle-Of-Attack,  a  [Degrees] 


Figure  22  Comparison  o  fCn^oYN  C-coefficient  of  the  Linear  Lateral/Directional  Characteristic 

Quartic  for  a  Preliminary  Configuration  of  the  X-31 A  Aircraft. 


NADC>90048-60 


Figure  23  CnpmN  Stability  Design  Guide  as  Suggested  by  Reference  (47) 
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As  illustrated  with  F-1 8  aircraft  data  in  Figure  24,  the  aerodynamic  rolling  and 
yawing  moments  of  many  of  today’s  fighter/attack  aircraft  are  nonlinear  with 
sideslip  angle  at  leu’ge  angles-of^attack. 


Sideslip  Angle,  p  ( Degrees ) 

Figure  24  Nondimensional  Rolling  and  Yawing  Moment  Coefficient  Versus  Sideslip  Angle  for  the  F-1 8 
Aircraft  at  30  Degrees  Angle-of-Attack  (Reference  (43)) 


The  analytical  technique  of  linearizing  Grand  Cn  at  zero  degrees  of  sideslip 
angle  in  these  cases  can  lead  to  erroneous  stability  estimates.  Given  this  fact 
along  with  the  fact  that  most  fighter/attack  aircraft  departure  susceptibility 
degrades  for  asymmetric  flight  conditions  makes  asymmetric  flight  an  impor¬ 
tant  consideration  in  the  determination  of  an  aircraft's  full  envelope  departure 
characteristics.  As  will  be  shown  later,  Pelikan’s  Cnp^pp  criteria  utilizes  a 
secant  slope  technique  (as  opposed  to  the  tangent  slope)  to  account  for  the  non¬ 
linearity  of  the  data.  STI,  Kalviste  and  Chody  on  the  other  hand,  linearize  the 
equations-of-motion  assuming  non-zero  sideslip  and  calculate  the  beta  deriva¬ 


tives  (C(p,  Cnp,  etc)  using  the  local  tangent  slope  (i.e.,  dCn/dp 


OO'XXX 

po=yyy 


). 
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(3)  The  absence  of  accounting  for  dynamic  maneuvering  effects 

n  “  1 1 

( i.  e.  ,p  =  Lp  +  -  +  •  ^Qor  +  RoqJ )  and  the  associated  longitudinal- 

lateral/directional  cross-coupling  as  disoissed  earlier  in  table  III. 

(4)  The  Cn^QYM  parameter  can  not  predict  a  divergence  that  is  caused  by  the  con¬ 
trols. 


At  high  angles-of-attack  it  is  very  common  for  today’s  fighter/attack  aircraft  for 
the  rudder  to  be  in  a  reduced  dynamic  pressure  area  of  the  wing  and 
forebody  wake  resulting  in  reduced  directional  control  effectiveness.  Wind 
tunnel  data  for  a  preliminary  X-31 A  configuration  examplifies  this  fact.  As 
shown  in  figure  25,  the  rudder  control  effectiveness  essentially  reduces  to 
zero  above  fifty  degrees  angle-of-attack. 


62 


NADC-90048-60 


At  the  same  time  rudder  control  effectiveness  is  degrading  with  increasing  angle- 
of-attack,  the  deflection  of  ailerons  at  high  angles-of-attack  often  results  in 
large  adverse  yawing  moments.  Using  the  same  X-31 A  configuration  as  an 
example,  figure  26  shows  that  differential  flap  deflection  at  low  angles-of-at¬ 
tack  produces  proverse  yaw,  whereas  at  the  higher  angles-of-attack,  above 
60  degrees,  full  differential  trailing-edge  flap  deflection  is  very  effective  in 
generating  adverse  yaw. 
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Figure  26  The  Effect  of  Differential  Trailing-edge  Rap  Deflection  on  the  Generation  of  Adverse  Yaw  at 
High  Angle-of  Attack  for  a  Preliminary  Configuration  of  the  X-31  A  Aircraft  (Reference  (11)). 


This  type  of  divergence  can  best  be  revealed  by  calculating  the  Lateral  Control 
Departure  Parameter  (LCDP)  discussed  next. 
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(5)  One  final  factor  pointed  out  in  Reference  (26)  which  can  affect  the  correlation 
of  the  parameter  with  actual  in-flight  departure  is  the  use  of  static  data 
which  do  not  accurately  match  flight  conditions  in  calculating  CnpoYn- 
most  influential  conditions  which  may  not  be  matched  are  control  deflections 
and  engine  power  setting.  In  some  cases  these  parameters  can  have  sub¬ 
stantial  effects  on  the  static  lateral/directional  stability  parameters 
Cnp  and  C|p. 
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BASIS  OFLCDP  CRITERIA: 

The  LCDP  parameter  is  derived  based  on  determining  conditions  for  which  the  roll  attitude  to 
aileron  control  input  transfer  function  ((ti(s)/6a(s))  is  nonminimum  phase  (i.e.,  when  the 

numerator  polynominal,  N$.,  contains  a  zero  in  the  right  half  plane  (RHP)  of  the  root  locus).  The 
equations-of-motion  utilized  are  the  linear  uncoupled  lateral/directional  dynamics  (i.e.,  rectilinesu’ 
flight,  <1)0-0  degrees  and  Vo  -  0  ft/sec)  assuming  the  reference  equilibrium  state  has  zero  an¬ 
gular  rates  (i.e.,  Po  -  Qo  -  Ro  -  0  deg/sec  and  vi»o  =  6o  =  ^  -  0  deg/sec). 

Routh's  "Stability"  Criteria  is  applied  to  the  polynominal  of  equation  (28)  to  determine  whether 
the  numerator  polynominal,  N^,  of  the  <Ks)/8a{s)  contains  any  RHP  zeros. 

N?,  •=  +  C4,  (28) 


The  numerator  polynominal  N$.  can  be  determined  by  taking  the  determinant  of  the  matrix 


given  by  equation  (29). 


S-A11  -Ai2 

-Ai3 

V6. 

-A21  S-A22 

-A23 

u. 

-A3I 

-A32 

S-A33 

Na. 

-A41 

-A42 

-A43 

0 

(29) 


Anm  -  The  elements  of  the  decoupled  lateral/directional  body-axis  Stability  Matrix  given  in  Figure  17. 


According  to  Routh's  Stability  Criterion,  the  quadratic  polynominal  given  by  N^,  is  guaranteed  to 
have  no  RHP  zeros  if  each  of  the  coefficients,  64.,  and  C^  as  given  in  equation  (28)  have 
positive  values.  (Note,  in  addition  to  the  assumptions  of  Figure  17,  the  algebraic  expression  for 
A<t„  B4,,  and  C4,  assumes  Yp  -  Yr  «■  Ya,  -  0). 


A4,  =  1 


+ 


•  tanOo 


(30a) 
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(30c) 


Moul  &  Paulson  found  that  negative  values  of  the  coefficient  were  associated  with 
lateral/directional  departures.  In  the  derivation  of  the  AADP,  the  condition  for  Cd,  having  a  posi¬ 
tive  value  was  simplified  to  the  expression  given  in  equation  (31  a)  by  assuming  that  theYv  term 
was  second  order  and  simplifying  the  first  term  by  assuming  the  use  of  stability  axis  (where 
U,;  =  Vt  :  Wo  =  0,  and  0o  =  to  which  is  assumed  to  be  zero.  Also  recall, 

N'p  =  N'v  •  Vt  :  af'p  =  a:'v  •  Vt)  yields. 


fN's 

Co,  =  (N'p)j-(a'p)^  ^  >0 

.  *  Js 


=  AADP>0 

or  expressed  in  nondimensional  terms, 

AADP  .  (C'np)^  -  (C'fp)^  7^  (31b) 

\  /$ 

Moul  and  Paulson  extended  the  concept  of  the  AADP  departure  susceptibility  parameter  to  also 
be  applicable  to  aircraft  with  an  aileron-to-rudder  interconnect  (ARI).  This  expression  is  given  by 
equation  (32). 
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DaForno  (Reference  (41 ))  offers  an  alternative  form  of  the  ARI-LCDP  parameter  that  allows  for 
a  more  direct  determination  of  the  &a/5r  crossfeed  necessary  to  maintain  LCDP  positive  (Refer¬ 
ence  (37)).  This  expression  is  given  by  equation  (33). 


•np 


r-  rr'-  _  n-  -v  n 

Ctj,  Cct,  Cjj,  6a 


(33) 


(Note,  all  derivatives  are  defined  in  the  stability-axis) 

Additionally,  Weissman,  in  Reference  (27)  provides  an  expression  for  LCDP  that  addresses  sys¬ 
tems  augmented  with  aileron  plus  rudder  proportional  to  sideslip  angle.  This  expression  is  given 
by  equation  (34). 


where;  Ki  •  - 


Cnj  C(p|^ 
6r 


rc, 


ns. 


+  Ki 


rc 


ns. 


7^Cc^-C 


ns, 


>  0 


(34) 


(Note,  all  derivatives  are  defined  in  the  stability-axis) 


TYPE.OE.CB1TER1A : 

Closed-loop  Lateral/Directional  Static  Stability  Criteria  with  respect  to  roll  control. 


CBlIEBlQIiPLOT: 


LCDP  DESIGN  IMPACT : 


The  LCDP  departure  susceptibility  parameter  primarily  indicates  conditions  under  which  roll 
reversal  (roll  is  in  the  opposite  direction  to  that  commanded)  is  likely  to  occur  and  divergence  or 
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departure  from  controlled  flight  is  likely  to  follow.  This  can  be  discerned  from  equation  (31).  The 
typical  condition  for  LCDP  having  a  negative  value  arises  when  yaw  due  to  "aileron"  (or  lateral 
control)  becomes  sufficiently  adverse  (negative)  that  the  second  term  overpowers  moderately 
positive  values  of  Cnp  (Reference  (37)). 

A  dynamic  system  becomes  unstable  in  a  closed-loop  sense  when  open-loop  system  poles 
(airframe  or  flight  control  system)  are  driven  toward  open-loop  zeros  located  in  the  RHP  (called 
nonminimum  phase  zeros)  of  the  o  -  jo)  root  locus  plane.  Nonminimum  phase  zeros  can  have  a 
significant  influence  on  an  aircraft's  transient  response.  Typical  of  nonminimum  phase  systems 
is  the  associated  control  reversal  of  the  initial  response  as  illustrated  in  figure  (27). 


Figure  27  Typical  Step  Responses  for  Nonminimum  Phase  Systems 
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The  nonminih.um-ph.se  ot  response  is  undesirable  (i.e.,  the  magnitude  «.d  sign  of  the  mo¬ 
tion  is  not  predlotabie)  because  olosed-loop  pilot  control  is  dilfloult  to  the  point  where  the  system 

can  be  driven  unstable. 

AS  mentioned  earlier,  aircraft  depahure  susceptibility  Is  both  an  open-loop  and  a  ctosed-loop 
phenomena.  As  witfi  most  flying  qualibes  criteria,  emphasis  ahould  be  placed  on  the  develop¬ 
ment  of  olosed-loop  parameters  and  relate  them  to  open-loop  paremmers.  In  the  study  ol  Beler- 
ence  (1 3)  one  ol  the  chlel  findings  was  that  departure  Is  -severely  aggravated  if  not  caused  by 
closed-loop  pilowehiole  interaction-.  Along  these  lines,  the  LCDP  is  the  most  rudimentary  and 
most  widely  accepted  closed-loop  departure  susceptibility  parameter. 
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WEISSMAN:  CnpovN  v®  LCDP 

DATE/RESEARCHER : 

1971;  WEISSMAN 

CB1I£BIQM : 

CnpoYN  versus  LCDP  Departure  Susceptibility  Criterion  Plane  as  shown  in  Figure  28. 
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Figure  28  Weissman  Copdyn  Versus  LCDP  Departure  Susceptibility  Criterion  Plane 


BASIS  OF  CRITERION: 

Weissman  developed  the  criterion  from  analyzing  time  history  sensitivity  studies  to  lateral/direc¬ 
tional  static  stability  derivatives  in  a  digital  six  degree-of-freedom  off-line  simulation.  Based  on 
these  time  history'  traces,  Weissman  empirically  identified  regions  of  increasing  roll  departure 
severity  and  spin  susceptibility  (see  Reference  (44)).  Weissman  later  correlated  flight  test  data 
(F-8,  F-102,  F-106  and  SAAB  37)  with  his  previously  defined  boundaries  (Reference  (45))  and 
found  good  agreement.  The  work  of  Titiriga  (Reference  (47))  modified  Weissman’s  original 
criterion  plane  by  adding  boundaries  (or  regions  E  and  F.  These  two  regions  delineate  suscep¬ 
tibility  of  an  aircraft  to  depart  controlled  flight  in  yaw  rather  than  roll  (see  Figure  29). 
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Figure  29  Weissman  LCDP  Versus  Departure  Susceptibility  Criterion  Plane  as  Modified  by  Skow 
and  Titiriga  of  Reference  (37) 

STI  has  in  turn  suggested  further  modifications  to  the  Northrop  Modified  Weissman  Departure 
and  Spin  Susceptibility  Criterion  Plane.  The  modifications  involved  shifting  the  E/F  boundary  to 
coincide  with  the  C/0  boundary;  raising  the  boundary  between  the  E/F  and  C/D  regions;  and  ex- 
tending  the  A/B  boundary  to  the  LCDP  axis.  These  changes  were  supported  by  the  piloted 
simulation  documented  in  Reference  (48).  The  STI  Modified  Weissman  Departure  and  Spin  Sus¬ 
ceptibility  Criterion  is  shown  in  Figure  30. 
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LCDP 
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Figure  30  Weissman  LCDP  Versus  CnooyM  Departure  Susceptibility  Criterion  Plane  as  Modified  by  STI 
(Reference  (48)) 


IYP£  OF  CRITERION: 

Combined  Open  and  Closed-loop  Lateral/Directional  Departure/Spin  Susceptibility  Criteria. 


CRITERION  PLANE : 

See  Figure  28,  29  and  30. 


DESIGN  IMPACT : 

The  authors  of  References  37  and  46  offer  an  insightful  design  interpretation  of  the  Cnpoy^  ver¬ 
sus  LCDP  criterion  plane  based  on  the  assumption  that  active  roll  control’  is  being  attempted  by 
the  pilot  (and/or  flight  control  system)  and  that  the  Cnpoy^  and  LCDP  departure  parameters  can  be 


’  Note  that  for  augmented  aircraft  the  use  of  any  feedback  to  rudder/aileron,  with  the  exceptior, 
of  roll  or  roll  rate,  will  alter  the  roll-  numerator  roots,  ui),,  and  hence  the  effective  LCDP  value 
(Reference  (37)). 
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approximated  by  the  simplified  Dutch-roll  pole-zero  frequency  relationships  of  equations  (35) 
and  (36). 


tOnoR ' 

llg  “•'Adyn 

(35) 

2-J 

(04,  *  - 

^LCDP 

(36) 

Figure  31  shows  the  pole-zero  relationships  associated  with  each  of  the  versus  LOOP 
criterion  plane  quadrants.  The  arrows  show  the  direction  of  the  closed-loop  root  migration  with 
increased  loop  gain. 

Region  A  depicts  three  cases  all  of  which  are  open-loop  stable  and  closed-loop  stable  at  high 
gain.  T  ne  differences  between  the  three  cases  is  as  follows: 

(1)  For  (04,/conoR  greater  than  one  the  locus  beginning  at  the  Dutch-roll  pole  will 
proceed  counterclockwise  into  the  zero  in  a  roughly  circular  arc.  Note  that  if 
^4,  and  Cor  are  lightly  damped  (i.e.,  dose  to  the  jcD-axis)  the  locus  may  pass  very 
dose,  or  even  into  the  unstable  RHP. 

(2)  For  (13,^/wnop,  equal  to  one  the  Dutch-roll  mode  is  not  excited  by  the  effective 
lateral  control  law  and  is  thus  decoupled. 

(3)  For  b}4,/(OnoR  less  than  one  the  locus  proceeds  counterclockwise  from  the 
Dutch-roll  pole  to  the  zero,  as  a  result  the  Dutch-roll  damping  increases 
slightly. 

Region  B  indicates  that  roll  control  reversal  exists  and  that  with  increasing  pilot  or  FCS  gain  the 
Dutch-roll  mode  will  eventually  go  unstable. 

Regions  C  and  D  indicate  an  airframe  instability  (Cnp^^^  <  0)  that  can  not  be  stabilized  by  roll 
control;  and  Regions  E  and  F  indicate  an  airframe  instability  that  can  be  compensated  via 
roll  control.  To  summarize,  when  LCDP  (ci>|)  is  negative,  closing  the  roll  loop  may  lead  to 
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Figure  31  Generalized  Pole-Zero  Relationships  Assodated  With  Quadrants  Of  The  LCDP  Versus 
^"foYN  Criterion  Plane  (Reference  (46)) 
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divergent  motions  or  if  Con>C(> «  ^  and  (04,/(onoR  >  1  closing  the  roll  loop  may  lead  to  oscillatory 
instabilities.  When  LCDP  is  negative,  the  possibility  of  stabilizing  the  Dutch-roll  mode  with 
aileron  control  vanishes  and  other  means  must  be  utilized  such  as  an  ay  or  p  feedback  to  rudder 
control  is  required. 

This  points  out  that,  when  designing  for  "good"  closed-loop  flying  qualities,  modifications  to 
basic  airframe  dynamics  as  well  as  the  PCS  design  to  obtain  favorable  zero  locations 

(i.  e.,  b>|  °  LCDP)  is  just  as  important,  if  not  more  important  than,  the  open-loop  eigenvalues. 


Prom  a  configuration  perspective,  low  or  negative  values  of  LCDP  can  be  improved  by  reducing 
the  effective  dihedral,  Cr^  and  or  utilizing  a  blended  spoiler/aileron  control  to  mitigate  adverse 
yaw  (Cnj,)  characteristics. 


Additionally,  via  flight  control  system  augmentation,  directional  stability  can  be  improved  by  aug¬ 
menting  Cnp  and  adverse  yaw  can  be  reduced  by  use  of  an  aileron-to-rudder  interconnect  con¬ 
trol  law  scheme.  Both  of  these  means  of  augmentation  require  the  availability  of  sufficient 
directional  control  power,  which  is  not  typically  available  at  high  angles-of-attack  using  conven¬ 
tional  aerodynamic  controllers.  The  X-31 A  aircraft  is  a  good  example  of  this.  As  illustrated  in  Pig- 
ure  25,  the  aerodynamic  rudder  effectiveness  is  effectively  zero  above  45  degrees 
angle-of-attack.  Above  45  degrees  angle-of-attack  the  yaw  thrust  vectoring  controller  must  be 
relied  upon  to  provide  all  the  directional  control  power  necessary  to  maintain  stability  and  control 
to  prevent  the  possibility  of  departing  controlled  flight. 


Besides  thrust  vectoring  control,  another  alternative  means  of  providing  directional  control  povrar 
that  is  beginning  to  show  some  promise  at  high  angles-of-attack  is  the  use  of  active  forebody 
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strakes.  The  strakes  are  used  to  control  the  relatively  povtrerful,  forebody,  vortices  that  eminate 
from  the  apex  of  the  sharp  forebodies  of  most  current  fighter  type  aircraft.  Reference  (49)  and 
(50)  are  good  sources  to  gain  more  information  concerning  current  research  advances  of  this 
new  configuration  design  technology. 

versus  LCDP  CRfTERiON  ADVANTAGES  ; 

The  use  of  the  Cn^Qy^  versus  LOOP  Departure  Susceptibility  Criterion  Plane  provides  two  major 
advantages  over  either  of  the  two  departure  parameters  or  LDCP  independently.  First  off 
it  provides  both  open  and  closed-loop  lateral/directional  stability  information  on  a  single  criterion 
plane  and  secondly  the  two  parameters  have  been  correlated  with  flight  test  data  to  identify 
regions  of  the  criterion  plane  in  terms  of  the  type  and  severity  of  departure  that  can  be  expected. 
A  third  and  important  advantage  of  this  criterion  is  that  the  patameters  can  be  calculated  early 
enough  in  the  design  phase  to  affect  the  final  configuration. 

SHORTCOMINGS  OF  THE  ^nppYN  versus  LCDP  CRITERION ; 

As  pointed  out  earlier  under  the  sections  that  discussed  the  Cppoy^  and  LCDP  departure 
parameters,  their  derivations  are  based  on  the  assumptions  of  uncoupled  longitudinal  and 
lateral/directional  dynamics  for  rectilinear,  steady  state  flight  conditions.  Additionally,  both 
parameters  neglect  second  order  terms  that  can  be  significant  for  certain  aircraft  configurations. 
This  section  briefly  summarized  the  evolution  of  the  Cnpjjy^,  versus  LCDP  Departure  Suscep¬ 
tibility  Criterion  and  one  should  observe  (as  pointed  out  by  MIL-STD-1797A)  that  using  the  same 
"static-based"  parameters,  independent  investigations  have  proposed  different  criterion  boun¬ 
daries  and  not  a  common  one.  This  points  out  that  one  must  clearly  understand  the  assump¬ 
tions  and  limitations  of  the  parameters  and  the  criterion  planes  based  on  these  parameters 
when  applying  them  to  "new"  aircraft  configurations  and  flight  conditions  that  violate  the  underly¬ 
ing  assumptions  and  data  base  used  to  develop  the  criterion. 
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JOHNSTON:  l/Teg 

DATE/RESEARCHEB: 

1974 :  Johnston  (STI ) 

CBITERION : 

No  departure  susceptibility  criterion  was  proposed  but  open-  and  closed-loop  departure 
parameters  were  identified  for  the  A-7  aircraft  considering  the  effects  of  nonzero  sideslip  angle 
and  the  coupled  ( longitudinal-lateral/directionsJ )  equations-of-motion.  The  closed  loop  nose- 
slice  departure  parameter  was  identified  as  1  /Tea  which  is  the  nonminimum  phase  zero  of  the 

N^M  transfer  function  of  the  complete  six  degree-of-freedom  bare  airframe  system  as  given  in 
equation  37. 


A 


Ke( -  1  /TeaX  1  /Te2)j^Cei.  ®ne,  JCe*.  conejJ 
4  Kp.  UnJ  ICrs*  “nnsl  Kdr>  Wnonl  Ksp>  ®nspl 


(37) 


The  nose  slice  departure  parameter  was  found  to  be  strongly  related  to  the  open-loop 
aerodynamic  parameters  N'a  and  ^'a  and  the  kinematic  terms  Zp(^  pocosao)  and  ZK^Posinao) 
as  given  in  Figure  32. 

In  trying  to  understand  and  predict  the  nose-slice  departure  characteristics  of  the  A-7  aircraft 
which  evolve  from  unsymmetrical  flight  (dire^ional  mistrim  or  miscoordination  during  maneuver¬ 
ing),  Johnston  concluded  that  although  pilot  ratings  (from  a  fixed  base  simulation)  could  not  be 
correlated  with  any  single  flying  quality  parameter,  the  A-7  departure  characteristics  were  most 
dependent  upon  the  coupled  longitudinal-lateral/directional  closed-loop  parameter,  1  /Tea  •  At 
the  A-7’s  departure  angle-of-  attack  the  parameter  1  /Tea  located  in  the  RHP  of  the  root  k>cus 
and  with  "normal"  pilot  pitch  attitude  control  activity  the  aircraft/pilot  system  is  driven  unstable. 
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Rgure32  Linear,  Rigid,  Six  Degree-of-Freedom  Coupled  ( 0  )Equations-of-Motion  for  the  Unuaugmented  A-7  Aircraft 
(Data  For  ao  -  18.8°;  (to  “  6°)  (Reference  (13)) 
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Figure  33  shows  a  comparison  between  the  bare  airframe  pole-zero  locations  of  the  un¬ 
coupled  and  coupled  N§,/a  transfer  function  for  the  A-7  aircraft  near  its  stall  angle-of-at- 
tack. 


:  6  DOF  coupled 


Figure  33  Comparison  Between  the  Pole-Zero  Locations  of  the  N&,M  Transfer  Function  for  the  Un¬ 
coupled  and  Coupled  6  DOF  Equation-of-Motion  for  the  A-7  Aircraft  ( ao  '  18.8®  ,  Po  “  6® ) 
(Reference  (13)) 


In  the  uncoupled  six  degree-of-freedom  case  the  two  lateral-directional  modes  ((i>non>  (^hsr)  have 
cancelling  pole-zero  dipoles;  the  longitudinal  zero  l/Tsj  is  small  and  positive  (LHP);  and  as 
might  be  expected  near  stall,  1  /Te,  is  slightly  negative  (RHP).  For  the  completely  coupled  six 
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degree-of-freedom  case,  the  coupling  has  a  relatively  small  effect  on  the  eigenvalues^  com¬ 
pared  to  a  major  shift  it  produces  on  the  zero  locations.  Of  primary  concern  is  the  shift  of  the 
1  /Te,  zero  into  the  RHP.  Also  of  concern  is  the  increased  separation  of  the  Dutch-roll  dipole 
pair  vtrhich  indicates  considerable  Dutch  roll  modal  response  excitation  exists  via  the  elevator 
controller. 

After  eliminating  the  U  equation  from  the  complete  equations-of-motion  and  observing  that  the 
zeros  identified  as  1  /Te,  and  1  /Tej  remained  unchanged  (the  1  /Te,  zero  is  usually 
eliminated)  in  location  from  the  complete  six  degree-of-freedom  case  (see  Figure  34),  Johnston 
deduced  that  the  nonminimum  phase  zero  is  a  result  of  the  lateral-longitudinal  coupling  and 
thereafter  identified  it  as  1  /Tea  (as  opposed  to  1  /Te, ). 

A  root  locus  showing  a  pure  gain  closure  of  pitch  attitude  to  elevator  is  shown  in  Figure  35  for 
the  six  degree-of-freedom  coupled  airframe  in  nonsymmetric  flight. 

Equation  38  provides  a  simplified  expression  for  1  /Tea  after  all  the  damping  derivatives  are 
neglected. 


1  /Tea  = 


Sf'gNp  -  Sf'pN'g 
sf'gcosao  +  N'gSinao 


(38) 


To  avoid  a  potential  closed-loop  "spiral"  divergence,  1  /Tea  '^ust  remain  positive. 

^  Nonzero  sideslip  angle  has  little  effect  on  the  longitudinal  eigenvalues.  The  major  shift  is  on 
the  lateral/directional  modes  in  which  increasing  sideslip  causes  the  roll  subsidence  and  spiral 
modes  to  couple  into  a  lateral  oscillation  (i.e.,  lateral  phugoid).  Typically  with  increased  sideslip 
angle,  the  Dutch-roll  and  lateral  phugoid  modes  interchange  damping.  That  is,  with  increasing 
sideslip  angle,  damping  of  the  Dutch-roll  mode  increases  while  the  damping  of  the  lateral 
phugoid  mode  decreases. 


80 


NADC-90048-60 


.0 


Nje  :  6  OOF  Coupled 


Figure  34  Comparison  between  the  Pole-Zero  Locations  of  the  N?,/4  Transfer  Function  for  the  6  DOF 
Equations-of-Motion  and  the  5  DOF  Equations-of-Motion  (No  0-equation)  for  the  A-7  Aircraft 
(  ao  “  18.8®  .  po  =  6®  )  (Reference  (13)) 
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Figure  35  Root  Locus  Plot  for  the  Six  Degree-of-Freedom  Coupled  0^0)^  Transfer  Function  (Data  for 
the  A-7  Aircraft;  ao  "  18.8®  ,  Po  *  6® )  (Reference  (13)) 
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The  root  locus  plot  shows  how  the  roots  starting  at  consn  rapidly  nnove  to  the  real-axis  and  then 
split  into  two  real  roots,  one  of  which  moves  towards  the  1  /Te2  zero;  the  other  moves  toward 
the  nonminimum  phase  zero,  1  /Tea . 


Figure  36  shows  the  system  time  response  when  the  pitch  attitude  loop  is  closed  with  unity  gain 
and  a  step  attitude  command,  Be,  of  0.01  radians  is  introduced.  The  aircraft  is  initially  trimmed 
for  steady  flight  at  oo  ■  18.8  degrees,  6o  ■  6  degrees,  and  410  ■  5  degrees.  The  predicted  first- 
order  divergence  is  shown  to  start  immediately  and  to  dominate  the  r,  v,  6  and  a  traces. 

The  cause  of  the  RHP  zero.  1  /Tea.  was  traced  back  through  the  9  x  9  equations  of  motion 
matrix  (see  Figure  32)  to  the  dominant  coupling  terms  ^'a,  N'a  ,  Zp  and  Zr.  ^'a  and  N'o  are  the 
dimensional  aerodynamic  cross-coupling  stability  derivatives  as  given  by  equations  39  and  40. 


a  = 


da 


gSb  .  dCt  .  aCn 

w,-& 


(39) 


N' 


a  ® 


iL 

da 


gSb 
IxU  “  1x2 


dCf  dCn 


(40) 


Wind  tunnel  tests  and  tuft  studies  performed  on  the  A-7  aircraft  revealed  that  the  decrease  in 
directional  stability  is  attributed  to  the  vortex  activity  from  the  fuselage  and  wing  center  panel  im¬ 
pinging  on  the  downward  side  of  the  vertical  tail  (Reference  (1 3)).  It  is  surmised  by  the  authors 
of  Reference  (13)  that  this  aerodynamic  phenomena  is  related  to  the  and  N'a  stability  deriva¬ 
tives  and  hence  1  /Tea-  The  effect  of  not  including  the  Jf'o  and  N’a  in  the  complete  linear  six  de- 
gree-of-freedom  equations  of  motion  of  the  A-7  can  be  inferred  from  Figure  37  which  compares 
the  pole-zero  locations  of  the  coupled  six  degree-of-freedom  N|./A  transfer  function  with  and 
without  the  and  N'a  terms. 
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Figure  36  Digital  Simulation  of  a  Closed-Loop  A-7  Departure  to  a  Step  Pitch  Attitude  Command 
(  ao  =  18.8®  ,  p  =  6®  )  (Reference  (13)) 
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N||,/A  ;  6  DOF  coupled 
(STa  -  N'a  -  0) 


Figure  37  tffect  i '  Not  Including  the  if  a  and  N'a  Terms  (in  the  6  DOF  Equations-of-Motion)  on  the  Pole- 
Zero  Locations  of  the  N?,/A  Transfer  Function  of  the  A-7  Aircraft  (Data  for 
«o  =  18.8®  :  Po  =  6®  )•  (Reference  (13)) 


The  Zp  and  Zr  terms  represent  the  nonlinear  kinematic  coupling  between  the  sideslip  and 
yawing/rolling  motion  of  the  aircraft  in  the  Z-equation-of-motion.  The  effect  of  not  including  these 
two  terms  in  the  six  degree-of-freedom  equations-of-motion  is  shown  in  the  pole-zero  plots  of  fig¬ 
ure  38. 
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N||,/A  :  6  DOF  coupled 
(Zp  -  Zf  -  0) 


Figure  38  Effect  of  Not  Including  the  Zp  and  Zr  Terms  (in  the  6  DOF  Equations-of*Motion)  on  the  Pole- 
Zero  Locations  of  the  6, t.  Transfer  Function  of  the  A-7  Aircraft  (Data  for  ao  '  18.8®  ;  Po  =  6® ). 
(Reference  (1 3)) 


TYPE  OF  CRITERION  : 

Values  for  the  closed-loop  parameter,  1  /Tea  open  loop  parameters  ^'a.  N'o,  Zp  and  Zr 

have  not  been  correlated  directly  with  departure  susceptibility  boundaries  to  be  considered  a 
"Criterion”.  However,  the  design  impact  and  influence  these  parameter  have  In  choosing  a 
methodology  to  predict  under  which  flight  conditions  an  aircraft  will  depart  controlled  flight  is  sub¬ 
stantial. 

CRiTERION  PLANE: 

The  nose  slice  departure  parameter,  1  /T©3  is  a  closed-loop  stability  parameter.  If  the  value  of 
1  /Tea  positive  a  roll-spiral  divergence  (as  described)  cannot  occur.  However,  positive  value 
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boundaries  have  not  been  correlated  with  aircraft  departure  susceptibility.  Conservatively  the 
value  of  1/Te3  should  be  positive. 

PESIQN  IMPACT : 

From  a  design  and  analysis  perspective  the  nose-slice  departure  parameter  points  out  that  there 
are  cases  in  which  aerodynamic  (longitudinal-lateral/directional)  coupling  becomes  a  significant 
influencing  factor  in  causing  an  aircraft  'o  depart  controlled  flight.  The  influence  of  these  terms  is 
not  discernable  in  departure  parameters/criteria  that  are  derived  based  on  uncoupling  the  lon¬ 
gitudinal  and  lateral/directional  equations  of  motion  (i.e.,  CnoDvrj  and  LCDP). 

In  fact,  for  the  A-7  aircraft,  the  authors  of  Reference  (1 3)  found  that  although  low  departure 
parameter  values  of  Cn^Qy^  and  LCDP  may  contribute  to  the  departure  of  the  A-7,  neither  is  a 
primary  factor.  However,  when  CnpovN  and  LCDP  were  negative  at  the  same  time  at  which  1  /Tos 
was  nonminimum  phase,  the  rapidity  and/or  severity  of  the  nose  slice  departure  is  increased 
(Reference  (13)). 

ADVANTAGES  OF  CRITERION : 

The  1  /Tes  closed-loop  stability  parameter  provides  stability  information  concerning  the  potential 
of  a  closed-loop  divergence  occurring  due  to  coupling  (both  aerodynamic  and  kinematic)  be¬ 
tween  the  longitudinal  and  lateral/directional  modes  of  motion.  This  is  the  first  stability  parameter 
presented  that  can  predict  the  existence  of  such  an  instability. 

Derivation  of  the  departure  parameter  1  /Tea  recognizes  the  importance  of  the  effects  of  sideslip 
(either  intentionally,  i.e.,  rudder  maneuvering  or  unintentionally,  i.e.,  adverse  aileron  yaw, 
mistrim,  etc.)  with  respect  to  aircraft  stability  and  control  at  high  angles-of-attack.  In  addition,  the 
influence  of  the  aerodynamic  cross-coupling  (Mp,  N’a,  ^’a)  and  the  dynamic  derivatives  (Cn„  Crp, 
etc.)  is  accounted  for  in  the  calculation  of  l/Tej. 
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CRITERION  SHORTCOMINGS ; 

The  1  /Tea  stability  parameter  was  derived  and  used  to  explain  the  departure  characteristics  of 
the  A-7  aircraft.  It  is  not  meant  to  be  used  as  a  departure  parameter  exclusive  of  other 
parameters  such  as  CnnQVN.  LOOP,  etc.  Rather  it  should  be  used  to  compliment  departure 
parameters  derived  based  solely  on  the  decoupled  equations-of-motion  to  determine  the  poten¬ 
tial  for  closed-loop  instabilities  arising  due  to  longitudinal-lateral/directional  dynamic  coupling. 
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KALVISTEi  Crriacop*  ^*^pcop>  ^ 

PATE/RESEARC.HEB : 

1 978;  Kalviste 

CRITERION  : 

Coupled  Parameters  Criteria, 

Npcop  ^  0 
MacoP  ^  0 
K<  1 

Kalviste's  coupled  parameters,  Cn^cop  •  Cmacop  ®''®  stability  parameters  that  indicate  the 
open-loop  static  rotational  stability  of  an  airaaft.  The  "cop"  subscript  stands  for  "coupled"  indicat¬ 
ing  that  the  parameters  are  derived  from  the  coupled  equations-of-motion  and  thus  are  ap¬ 
plicable  to  asymmetric  flight  conditions. 

BASIS  OF  CRITERION; 

Kalviste's  stability  parameters  are  derived  based  on  solving  for  the  "necessary"  conditions  for 
aircraft  static  rotational  stability  (note,  the  aircraft  may  still  be  dynamically  unstable)  for  the 
coupled  six  degree-of-freedom  equations  of  motions  by  applying  the  Routh  Stability  Criterion.  A 
detailed  development  of  Kalviste's  parameters  from  the  equations-of-motion  is  given  in  Refer¬ 
ence  (51)  and  will  not  be  repeated  here.  Instead  the  implications  and  understanding  of  these 
parameters  will  be  elaborated  on. 

The  nondimensional  coupled  stability  parameters  are  defined  in  equations  (41)  to  (47). 


Npcop  ”  2 

1 

^NpoYN  ~  ^NpoYN  '*■  ^ODYNj  ■  i 

*  2 

\ 

1  ^NpoYN  +  MqoyN  j 

(41) 
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M. 


acop 


^NpDYN  ~  ModYN  j  +  MqdyN  j 


^  2 
^NpDYN  ■*■  MoOYNj 


(42) 


..A 

K  -  2 

^NpoYN  ■*■  ^^aOYNj 


(43) 


where, 


MaoYN  “  Mo  -  (^lacosao  +  NaSinao)tanPo 


(44) 


MpoYN  =  Mp  -  (Sfpcosao  +  Npsinoojtanpo 


Nodyn  =  No  cosao-  SfoSinao 


NpoYN  “  Np  cosao”  Sfpsinao 


(45) 

(46) 

(47) 


The  equations-of-motion  used  in  the  derivation  of  the  Kalviste  stability  parameters  differ  in  two 
major  ways  from  those  used  to  derive  CnpoYw  • 

1 .  The  aerodynamic  cross-coupling  between  the  longitudinal  and  lateral/direc¬ 
tional  equations-of-motion  are  modelled  (i.e.,  ^a.  Nq,  Mp) 

2.  The  kinematic  cross-coupling  terms  are  modelled  by  the 
(cosaotanpo)  and  (sinaotanpo)  terms. 

If  the  sideslip  angle  is  zero  then  Mqoyn  reduces  to  Mo  •  Furthermore  if  the  aerodynamic  cross¬ 
coupling  terms  are  zero  then  the  parameters  Nooyn  Mpdyn  are  both  zero  and  the  Kalviste 
Stability  Parameters  simplify  to  Npqyn  *  Np^op  >  0  and  Mo  -  Mocop  <  0- 
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Analogous  to  Npqyn  and  Mqoyn  being  equal  to  the  square  of  the  Dutch-roll  and  short  period 
natural  frequencies  respectively  (assuming  zero  damping;  see  equation  (48),  the  coupled  ex¬ 
pression  for  these  two  parameters,  NpcxM>  and  Mqcop  •  can  be  equated  to  the  coupled  Dutch-roll 
and  short  period  natural  frequencies  (see  equation  49). 

UNCOUPLED  (Mpov^NaovN  '  0)  (48) 

2 

^PoYN  “ 

^  2 

^aovKi  “  Wnsp 

COUPLED  ^MpoYN^aovN  *  Oj  (49) 

•  2 

Npoop  “  (“non)cop 

^Ocop  “  (wnsp)cop 


To  predict  aircraft  departure  susceptibility  for  a  particular  aircraft,  Kalviste  makes  use  of  contour 
mapping  techniques  as  shown  in  Figure  39  to  define  three  regions  of  unaugmented  aircraft  in¬ 
stability  as  a  function  of  angle-of-attack  and  sideslip  angle.  These  three  regions  of  instability  are; 

1 .  K  <  1  Coupled  (a,  P)  Oscillatory  Instability 

2.  Macop  >  0  Coupled  Longitudinal  Divergence  (i.e..  Short  Period  on  RHP  o-axis) 

3.  NpcoP  <  0  Coupled  Lateral- Directional  Divergence  (i.e..  Dutch-roll  root  on  RHP  o-axis) 


90 


NADC-90048-60 


et 

DEG 


pDEG 


rrn  k<o. coupled (b.P i 
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RB  Cnp  COP*  0  .  COUPLED  LATERAL/ 
DIRECTIONAL  DIVERGENCE 


5 


Figure  39  Kalviste  a  versus  p  Static  Stability  Plot  (Reference  (51)) 


If  the  angle-of-attack  and  sideslip  traces  of  an  aircraft  in  maneuvering  flight  pass  through 
regions  of  instability  it  indicates  that  the  aircraft  will  have  a  tendency  to  depart;  it  does  not  neces¬ 
sarily  mean  that  the  aircraft  Ml  depart  controlled  flight.  There  are  two  possible  effects  of  an  un¬ 
stable  region  on  aircraft  motion.  They  are:  (1)  if  the  unstable  region  is  small  it  can  diverge  into  a 
stable  region;  (2)  if  the  unstable  region  is  large,  the  divergence  can  cause  the  aircraft  rates  to 
build-up  into  a  developed  post-stall  gyration  or  spin. 


It  is  important  to  realize  that  regions  of  the  "Kalviste  stability  plane"  that  are  not  identified  as 
regions  of  instability  do  not  imply  that  departure  cannot  occur  in  these  regions.  Aircraft  departure 
from  controlled  flight  due  to  dynamic  instability  (i.e.,  as  a  result  of  inertial  coupling  effects,  or 
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negative  roll  and/or  yaw  damping)  or  closed-loop  instabilities  due  to  pilot  or  flight  control  system 
loop  closures  are  not  addressed  by  the  Kalviste  open-loop  Coupled  Stability  Parameters. 

TYPE  OF  CRITERION: 

Open-loop  longitudinal-lateral/directional  static  stability  criterion. 

gBlTERlQN  PLANE : 

See  a  -  p  Stability  Plane  of  Figure  38. 

DESIGN  IMPACT : 

Kalviste's  stability  parameter  analysis  technique  can  be  pursued  after  nonlinear  static  wind  tun¬ 
nel  data  becomes  available.  To  make  the  analysis  feasible  requires  the  use  of  digital  computer 
programming  and  automatic  plotting  capability.  The  computer  is  used  to  perform  nonlinear  inter¬ 
polations  of  tabular  function  (i.  e.,  Cf  -  /((a,  p,  6a, ),  Cm  ■/  ((M,  a,  p,  5c))  that  have  continuous 
first  derivatives  through  the  a,  p  range  desired.  An  iteration  procedure  is  then  used  to  compute 
the  stability  contour  lines. 

CRITERION  ADVANTAGES : 

The  advantage  of  Kalviste's  Coupled  Static  Stability  Parameters  criterion  is  that  it  extends  the 
open-loop,  symmetric,  static  Cnpoy^  stability  parameter  to  asymmetric  flight  (which,  as  shown, 
can  be  a  prime  driver  in  initiating  an  aircraft  departure)  and  includes  the  effects  of  highly  non¬ 
linear  aerodynamic  data  typical  of  high  angle-of-attack  flight  conditions. 

CRITERION  SHORTCOMINGS : 

The  Coupled  Static  Stability  Parameters  criterion  does  not  address  the  potential  for  departure 
due  to  dynamic  instabilities  (i.e.,  aerodynamic,  inertial  coupling,  etc.),  nor  does  it  address 
aspects  of  potential  departures  due  to  closed-loop  pilot  control. 
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DATE/RESEARCHER  : 

1978;  Bihrle 

CRITERION  : 

C(p  versus  Cno  and  Cnt,  Design  Boundaries  for  Departure  Susceptibility  and  Roll  Reversal  (see 
Figure  40) 

BASIS  QF-CRUEBLQN: 

Experience  led  Bihrle  to  believe  that  there  are  three  aircraft  characteristics  primarily  responsible 
for  departure  susceptibility.  These  characteristics  are  static  directional  stability  (Cn|i),  dihedral  ef¬ 
fect  (C(^)  and  lateral  control  Cn^,).  To  determine  boundaries  for  aircraft  departure  suscep¬ 
tibility  and  roll  reversal,  Bihrle  conducted  a  digital  simulation  of  a  severe  open-loop  rolling  pull-up 
maneuver  while  parametrically  varying  the  Cn^  and  aerodynamic  parameters  for  three  cases 

of  yaw  due  to  aileron  control  (i.e.,  proverse,  neutral  and  adverse).  The  digital  simulation  utilized 
the  complete  6  degree-of-freedom  rigid  body  equations-of-motion  incorporating  the  use  of  non¬ 
linear  aerodynamics  with  angle-of-attack  and,  if  appropriate,  control  deflection  and/or  sideslip 
angle.  The  choice  of  the  severe  rolling  pull-up  maneuver  for  the  digital  simulation  is  based  on 
the  fact  that  most  fighter  aircraft  departures  from  controlled  flight  occur  while  maneuvering.  The 
particular  maneuver  chosen  is  also  severe  in  terms  of  the  amount  of  kinematic  and  inertial  cou¬ 
pling  generated. 

By  examining  the  resulting  time  history  associated  with  each  of  the  matrix  variations  of 
Cnp  and  C(p  the  boundaries  shown  in  Figure  40  were  empirically  determined.  Bihrle  defined  roll 
reversal  as  roll  opposite  in  direction  to  what  was  commanded  by  lateral  control  input.  Departure 
was  considered  to  have  occurred  when  the  angle-of-attack  was  sustained  above  the  maximum 
trim  angle-of-attack  value. 
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Figure  40  Bihrle  Applied  Research  Design  Guidelines  for  Departure  Susceptibility  and  Roll  Reversal  ( Reference  ( 53) ) 
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Some  of  the  noteworthy  aesumptions  that  influenced  the  results  of  the  design  charts  of  Figure 
40  are  as  follows: 

1 .  The  parametric  variation  of  Cn^  and  Cr^  (nonlinear  with  angle*of*attack)  was  based  on 
data  representative  of  current  high  performance  fighter/attack  aircraft  at  the  time,  such 
as  the  F-4.  F-1 1 1 ,  F-14  and  F-15. 

2.  The  ratio  (Iz/Ix)b  was  fixed  at  a  value  of  6.2  and  ( I  xz)b  was  assumed  to  be  zero. 

3.  The  boundaries  were  developed  for  longitudinally  stable  aircraft  only. 

4.  There  was  no  limitation  placed  on  the  longitudinal  control  authority. 

IYRE.OF-CRJTEBIQN : 

Empirical  (based  on  current  high  performance  fighter/attack  aircraft  data  of  the  time)  static 
stability  design  charts  applicable  to  large  angle  transient  maneuvering  flight. 

CRITERION  PLANE ; 

See  Figure  40. 

DESIGN  IMPACT : 

With  respect  to  departure  resistant  aircraft  design  and  positive  roll  control,  Bihrle's  empirical 
based  "design  charts"  provide  similar  information  to  the  Weissman  Criteria  in  terms  of  desired 
aircraft  static  stability  characteristics.  That  is,  the  "ideal"  aircraft  static  characteristics  are  a  com¬ 
bination  of  positive  effective  dihedral  (negative  Ctp)  along  with  positive  static  directional  stability 
(positive  Cnp)  and  low  adverse  yaw  lateral  control  characteristics.  Bihrle's  design  charts  show 
that  smaller  values  of  positive  effective  dihedral  can  be  offset  by  increased  directional  stability. 
On  the  other  hand,  static  directionally  unstable  aircraft  configurations  can  still  be  made  depar¬ 
ture  resistant  if  the  effective  dihedral  is  sufficiently  large.  However,  for  this  situation  the  occur¬ 
rence  of  roll  reversal  is  more  likely. 
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CRITERION  ADVANTAGES  : 

The  principal  advantage  that  the  Bihrle  Departure  and  Roll  Reversal  Design  Charts  have  to  offer 
is  that  they  are  based  on  maneuvering  flight  conditions  and  thus  address  maneuvering  stability. 
By  proceeding  in  this  manner  the  effects  of  kinematic  and  inertial  coupling,  as  well  as  the  effects 
of  nonlinear  aerodynamics,  are  addressed. 

CRITERION  SHORTCOMINGS : 

The  primary  shortcoming  of  the  Departure  and  Roll  Reversal  Boundaries  developed  by  Bihrle  is 
that  they  are  both  a  function  of  the  open-loop  large  angle  nonlinear  maneuver  simulated  and  the 
aerodynamic  data  base  and  aircraft  configuration  (i.e.,  (1z/Ix)b  *  6.2,  etc.)  modelled.  The  author's 
of  Reference  (37)  comment  that  they  believe  the  Bihrle  Departure  and  Roll  Reversal  Boundaries 
are  conservative  based  on  the  severity  of  the  dynamics  associated  with  the  rolling  pull-up 
maneuver  simulated.  The  Bihrle  design  boundaries  have  been  correlated  with  the  A/D  sectors  of 
the  empirically  derived  Weissman  Criterion  Plane  in  Reference  (37).  However  Reference  (37) 
cautions  against  comparing  the  two  criteria  because  the  Bihrle  boundaries  are  based  on  fixed 
ratios  of  inertias  (Iz/lx)  and  aileron  moments  (Cfe,/Cn»,),  where  as  the  Cn^oYN  LCDP 
parameters  that  comprise  the  Weissman  Criteria  are  not  restricted  in  this  manner. 
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CRITERION : 

1  /T(^,  or  C<t>o>n«  >  ~  0.5  for  Departure  Resistartce  where  1  /T^,  is  the  non-minimum  phase  zero  of 
the  (3  DOF)  transfer  function  typical  of  fighter  aircraft  at  high  angle-of-attack  flight  condi¬ 
tions. 


When  the  conditions  are  satisfied  that  permit  the  longitudinal  and  lateral/directional  equations  of 
motion  to  be  decoupled,  the  3  degree-of-freedom  transfer  function  can  be  represented 


by  equation  50  given  below. 


N^stk  _ 1^  tk[s^  ♦  2C(»MrvS  ii)^| _ 

^  (S'*’  1/Ts)(S  1/Tb)^S^  +  ^CoR^rion  *  "*■  wnonj 


(50) 


Typically  for  low  angle-of-attack  flight  conditions  the  numerator  zeros  are  a  complex  pair  located 
in  the  neighborhood  of  the  Dutch-roll  poles  (i.e.,  °  °  unon)  bs  shown  in 

Figure  41. 


Figure  41  Typical  Low  Angle-of-Attack  <t)(s)/5sTK(s)  Root  Loci  (Reference  (1 3)). 


97 


NADC-90048-60 


However,  with  increasing  angle-of-attackthe  ton*  or  C-coefficient  ^N^stk  *  As^  +  Bs  +  Cj 

of  the  N^stk  numerator  polynominal  decrease'',  and  can  bec'^me  negative.  When  th.s  occurs  the 
numerator  polynominal  can  be  expressed  as  given  in  equation  51  where  one  root  is  positive  (by 
definition  the  1  /T,^,  root)  and  one  root  is  negative. 

^5stk  “  (51) 

In  this  case  the  root  locus  of  Figure  42  shows  that  a  roll  command  loop  closure  will  diive  the 
spiral  root  (-1  /T$)  toward  the  RHP  (nonminimum  phase  zero)  and  a  first  order  divergence 
results.  (Note:  the  shift  in  the  open-loop  denominator  roots  is  common  at  high  angles-of-attack). 


Figure  42  Typical  High  Angle-of-Attack  4(8)/&stk(s)  Root  Lod  (Referenoe  (13)) 

The  rate  of  divergence  depends  upon  how  far  the  zero  lies  in  the  RHP  and  how  tightly  the  roll 
loop  is  closed  (i.e.,  pilot  technique,  aggressiveness,  etc.).  In  Reference  (1 3)  the  closed-loop 
parameters  un,  and  (2Co,(on,)  or  (1  /T,>„  1  /T^,,)  were  plotted  versus  angle-of-attack  for  the  F-4J 
aircraft  as  shown  in  Figure  43. 


98 


NADC-90048-60 


Figure  43  N^stk  Root  Migration  with  Angle-of-Attackforthe  F-4J  Aircraft  (Reference  (13)). 
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This  figure  illustrates  two  important  points.  First,  note  (as  mentioned  above)  that  negative  values 
of  (which  the  roll  reversal  parameter  LOOP  is  an  approximation  of)  implies  the  numerator 
zeros  can  now  be  represented  by  two  real  roots,  one  of  which  is  nonminimum  phase.  The  first 
point  the  authors  of  Reference  1 3  make,  is  that  for  the  F-4J  aircraft  (at  zero  sidesiip  angle)  wing 
rock  and  eventual  nose  slice  characteristics  are  correlated  with  small/negative  values  of 
and  nose  slice  correlates  with  large  negative  values  of  (1  /T^,,). 


The  second  point  illustrated  is  the  effect  of  nonzero  sideslip  angle  (i.e.,  5.5°)  on  these 
parameters.  For  this  F-4  configuration,  aerodynamic  cross-coupling  causes  (i.e.,  LCDP)  to 
remain  positive  and  close  to  the  Dutch-roll  mode.  However  the  damping  ratio  becomes 
negative  causing  two  nonminimum  phase  zeros  to  exist  in  the  form  of  a  complex  pair  as 
sketched  in  Figure  43.  This  points  out  one  case  in  which  LCDP  fails  to  predict  the  existence  of  a 
nonminimum  phase  zero  due  to  the  affects  of  significant  aerodynamic  cross-coupling.  In  this 
case  the  parameters  1  /T^.,  or  correlate  more  closely  to  the  departure  characteristics  of  the 
F-4J  aircraft. 


In  the  piloted  simulation  reported  in  Reference  (48),  four  configurations  were  evaluated  in  terms 
of  their  departure  susceptibility  in  an  attempt  to  isolate  the  influence  of  the  open-loop  dynamics 
only  (i.e.,  Cppoyn)-  closed-loop  dynamics  only  (i.e.,  LCDP,  1  /T^„  Cion*)  and  the  combination  of 
degraded  open  and  closed-loop  dynamics.  The  results  of  the  simulations  as  supported  by  the 
data  of  Figure  44  shows  the  division  between  "departure  resistant"  (R),  and  "departure  suscep¬ 
tible"  (S)  ratings  to  lie  at  an  approximate  1/T^,  value  of  >0.5. 
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KEY:  O  pact  JF 
O  PBot  RC 


(1/sec) 

Closed  Loop  Divergence  Potential 


Figure  44  Departure  Susceptibility  Ratings  Versus  Lateral  Closed-Loop  Divergence  Potential, 

1  /T^,,  or  5^,  (On*  (R  -  Departure  Resistant;  S  -  Departure  Susoeptible;  ES  -  Extremely 

Departure  Susoeptible)  (Reference  (48)). 
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This  corresponds  to  a  time  to  double  amplitude  of  approximately  1.4  seconds.  The  authors 
reasoned  that  zeros  which  lie  to  the  left  of  this  line  apparently  limit  the  first  order  divergence  to  a 
rate  slow  enough  for  pilots  to  respond  and  prevent  departure  from  controlled  flight.  Zeros  to  the 
right  of  the  -0.5  boundary  apparently  allow  divergence  rates  so  fast  that  the  pilots  can  not 
prevent  departure.  Interestingly,  for  the  flight  conditions,  inertias,  etc.,  employed  in  this  simula¬ 
tion  a  1  /T0,  value  of  -0.5  corresponds  to  an  LOOP  value  of  -0.001 .  This  coincides  with 
Weissman's  LCDP  boundary  between  regions  A  and  B  for  positive  CnpQy^  (see  Figure  28),  how¬ 
ever  it  is  a  little  more  conservative  at  negative  Cnpoy,^  values  (Reference  (1 3)). 

TYPE  QFCBITERIQN: 

The  1  /Tj,,  criterion  is  intended  as  a  closed-loop  departure  parameter  in  terms  of  preventing  un¬ 
commanded  motion  via  roll  control  inputs.  Unlike  the  LOOP  parameter,  1  /T^,  is  directly  ap¬ 
plicable  to  augmented  aircraft. 

CRITERIQliELANE ; 

1/T<„>  -0.5 

DESISNIMEACJ: 

The  design  implication  is  that  if  the  combined  aerodynamics  and  FCS  design  is  such  that  the 
value  of  1  /T^,  (and/or  )  does  not  exceed  -0.5  (negatively)  throughout  the  achievable  angle- 
of-attack  and  sideslip  range,  the  aircraft  is  predicted  to  be  free  of  roll  reversal  characteristics. 

The  1  /T,>,  departure  parameter  addresses  the  influence  of  the  FCS  on  departure  and  can  be 
used  to  aid  in  the  design  of  departure  prevention  flight  control  system  architecture.  This  is  par¬ 
ticularly  important  for  departure  design  of  future  fighter  aircraft  which  will  very  likely  have 
reduced  open-loop  directionally  stability  (i.e.,  vertical  tailless,  low  observable  configuration 
designs).  Note  that  this  criterion  places  no  restrictions  on  open-loop  {i,e.,  Cnpoy^)  stability. 
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CRITERION  ADVANTAGES: 

Unlike  application  of  the  closed-loop  departure  parameter,  LCDP,  the  1  /T^,  parameter  is  not 
restricted  to  the  three  degree-of-freedom  iateral/directional  equation-of-motion  assumptions. 

The  parameter  can  be  applied  to  completely  coupled  six  degree-of-freedom  airframe  dynamics 
with  FCS  augmentation. 

CRIIEBIQRSHPRISOMmSS : 

Admittedly  the  researchers  make  clear  that  the  simulation  results  used  to  develop  the  1  /T(|„ 
criterion  are  based  on  variations  of  a  single  nonlinear  aerodynamic  model  (F-4J)  representing 
an  a,  p  region  dominated  by  phenomena  that  are  highly  configuration  dependent.  Therefore,  fur¬ 
ther  substantiation  of  this  criterion  is  desired. 
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PELIKAN:  Cnp^pp 

DATE/RESEARCHER ; 

1984;  Pelikan  (McDonnell  Douglas) 


CRITERION : 

Cnp^p(i.e.,  CnpDYN  .  CONTROLS^  ^  ® 

( ACn).  +  ( ACn).  (AC() .  +  (AC().  . , 

Where,  Cnj^pp  - - - cosa - ^-r— - sinal 


(51) 


This  parameter  has  been  used  to  assess  the  effect  of  different  control  inputs  (i.e.,  neutral  con¬ 
trols,  lateral  stick,  coordinated  controls,  cross  controls,  rudder  alone)  on  the  basic  airframe  static 
lateral/directional  stability. 


BASIS  OF  CRITERION; 

The  CnpApp  departure  parameter  is  based  on  whether  the  initial  stability-axis  yawing  moment  ac¬ 
ting  on  the  aircraft  due  to  both  sideslip  and  control  inputs  is  either  restoring  or  propelling.  The  ex¬ 
pression  for  CnpAPp  9'ven  in  equation  (51)  can  be  derived  by  transforming  the  body  angular 
accelerations  (N/Izg  yi/lxe)  from  either  body  or  principal  axis  into  stability  axis  using  a  secant 
slope  linearization  technique.  One  development  of  the  equation  that  defines  the  Cnp^p 
parameter  is  given  in  equations  (52)  to  (54). 

Rs  *  Recosa  -  Pesina  (52a) 

■  ■  1^'"®  (®2b) 


(52c) 


t^Sb 


•  Rj  =  CnCOSa 


CfSina 


(52d) 
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Defining  the  nondimensional  body-axis  yawing  and  rolling  moment  coefficients  with  regard  to 
the  static  terms  due  to  sideslip  (p)and  control  deflections  (6)  yields, 


■  Rs  '  I(ACn)p 


+  (ACn)j]  •  COSO  - 


•  I(ACf)p  +  (ACt)^]-  sino 


(53) 


To  arrive  at  Pelikan’s  final  expression  for  equation  (53)  is  divided  through  by  the  sideslip 
angle, 


''"Papp 


p(ACn)j  +  (ACn)^- 

rn-  fin 

r[qsbj 

[  P  J 

“  [si  P  J 

■  sino 


(54) 


Pelikan  s  motivation  for  dividing  equation  (54)  by  the  sideslip  angle  is  to  use  a  secant  slope 
linearization  method  (as  opposed  to  the  local  tangent  slope  linearization  method,  i.e.,  stability 
derivatives)  to  determine  "static  stability"  in  a  "global"  sense  as  a  function  of  sideslip  angle  (see 
Figure  45). 

Note  that  the  Cnp^pp  parameter  derivation  is  based  on  the  same  concept  as  Moul  and  Paulson's 
CnpuYN  parameter.  That  is,  the  summation  of  yawing  moments  about  the  Z-stability  axis  must  be 
positive  (i.e.,  ^  N's,  >0)  (Reference  (37)).  The  difference  lies  in  the  inclusion  of  yawing  and  rolling 

i 

moments  due  to  nonzero  control  deflections  and  the  use  of  a  secant  slope  linearization  technique. 
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dC' 

Tangent  Slope  (M  =  ,  i  •  C,  n) : 

ap 


S«antStop.(M-.gg^):. 


c/ 


Figure  45  Illustration  of  the  Tangent  and  Secant  Slope  Linearization  Techniques 


From  F/A-1 8  flight  test  data  Pelikan  observed  open-loop  instabilities  predicted  by  the 
departure  susceptibility  parameter  (Coppy^  <  0)  that  were  never  realized  in  flight.  Based  on  this 
observation,  Pelikan  felt  that  the  association  of  local  open-loop  instabilities  (small/negative 
values  of  Cnppy^ )  with  departure  susceptibility  should  be  revised  to  use  the  secant  slope  method 
of  calculatingCnppy^  when  nonlinear  lateral-directional  characteristics  are  exhibited.  In  the  F-18 
case  the  local  region  of  instability  occurring  at  approximately  30  degrees  tends  to  slide  the 
aircraft  out  to  a  stable  region  between  6  and  8  degrees  sideslip  angle  and  prevents  the  aircraft 
from  departing  controlled  flight  (although  diret^ional  controllability  is  sloppy  in  this  region). 
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TYPE  OF  CRITERION: 

Open-loop  lateral/directional  static  departure  susceptibility  criterion  that  includes  the  static  effect 
of  the  controls. 

CRITERION  PLANE: 

DESIGN  IMPACT: 

The  intent  of  the  Cnp^p  criterion  is  to  extend  the  application  of  the  Cnpoy^  parameter  to  include 
nonlinear  aerodynamic  characteristics  with  sideslip  angle  and  the  effect  of  different  control  in¬ 
puts.  Use  of  this  criterion  is  useful  in  preliminary  control  law  design  architecture  for  departure 
prevent  portions  of  the  flight  control  laws. 

CRITERION  ADVANTAGES: 

(See  Design  impact) 

CRITERION  SHORTCOMINGS: 

The  CnpApp  parameter/criterion  does  not  address  the  potential  for  departures  due  to  dynamic  in¬ 
stabilities,  (i.e.,  aerodynamic,  inertial  coupling,  etc.)  nor  does  it  address  aspects  of  potential 
departures  due  to  closed-loop  pilot  control.  Additionally,  the  parameter  exclusively  considers 
departure  only  in  the  lateral  directional  axis,  precluding  the  possibility  of  any  longitudinal  axis 
departures  or  longitudinal-lateral/directional  coupling  that  might  induce  a  departure.  It  is  also 
believed  that  further  validatiort/correlation  with  flight  test  data  is  needed  to  support  use  of  the 
secant  slope  method  of  linearization  that  is  the  basis  of  the  Cnp^  parameter  derivation. 
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KALVISTE:  DYNAMIC  STABILITY  PARAMETERS 
DATE/RESEARCHER : 

1989;  Kalviste  (Northrop) 

CBIJEBIQN: 

No  departure  susceptibility  criterion  was  proposed.  Rather,  *t)ew"  stability  parameters  are 
derived  (given  in  equation  (55)  below)  whidt  relate  directly  to  the  lasf  aircraft  stability  modes 
(pitch  short  period,  Dutch-roll  and  roll  subsidence). 

^^Mcop  ^  0,  (55) 

Npcop^  0 
NpcoP>  0 
MacOP  ^  ^ 

MocoP  ^  ® 

Failure  to  satisfy  the  stability  conditions  of  equation  (55)  reveals  regions  of  local  instability  which 
result  in  degraded  unaugmented  flying  qualities  such  as  unpredictable  control  response  and  the 
inability  to  transition  to  and  maintain  certain  Uim  conditions  (Reference  (5)}.  The  local  in¬ 
stabilities  however  may  not  necessarily  represent  conditions  vrhere  the  aircraft  will  depart  from 
controlled  flight. 

These  five  parameters  are  related  directly  to  the  three  (roll  subsidence.  Dutch-roll,  and  pitch 
short  period)  unaugmented  stability  modes  (the  six  eigenvalues)  of  the  aircraft  in  terms  of  three 
new  independent  variables  ( pcop,  Pqop  )  >>  •  function  of  p,  P  and  a.  As  yet,  Kalviste 

has  not  developed  a  unique  transformation  between  the  uncoupled  variables  (p,  fi  and  a)  and 

the  coupled  variables  ( pcop,  0oop  >nd  )  for  the  general  case  involving  both  static  and 
dynamic  stability.  A  transformation  between  the  uncoupled  variables  (P  and  a)  and  the  coupled 
variables  (  and  )  has  been  derived  in  Reference  (4)  for  the  reduced  problem  involving 
static  stability  only. 
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Th«  necessary  condition  for  aerodynamic  stability  is  to  have  the  roots  of  the  coupled  equations 
in  the  left  half  of  the  complex  plane. 

BASIS  OF  CRITERION: 

The  five  stability  parameters  given  in  equation  (55)  define  aircraft  stability  based  on  the  aircraft's 
aerodynamic  and  inertial  properties  and  include  both  static  and  dynamic  aerodynamic  effects,  in* 
ertial  coupling  and  kinematic  coupling  effects.  These  parameters  are  an  extension  of  the  pre¬ 
viously  developed  Coupled  Stability  Parameters  based  only  on  static  aerodynamic  effects  (see 
Reference  (48)).  Kalviste  defines  stability  in  terms  of  the  rotational  motion  relative  to  the  flight 
path.  The  aerodynamic  moments  and  forces  are  functions  of  the  aerodynamic  angles,  a  and  p, 

aircraft  body  axis  rates,  P.  O,  and  R,  and  translational  acceleration  terms,  a  and  0,  at  a  fixed 
flight  condition  (Mach,  altitude).  To  derive  expressions  for  the  stability  parameters  of  equation 
55,  the  rotational  equations  of  motion  are  first  expressed  in  the  dynamic  stability  axis  system 
(i.e.,  the  coordinate  system  about  which  a,  p,  and  p  are  defined).  These  are  shown  in  equations 
56to67^ 

p  “  *  Sf^OYNp  *  ^PoynP  *  SfeovN®  ■*“  S^oovnO 

P  ■  ~  NpQyuP  "  NpDYNP  “  NPqymP  -  NaoY|i|6l  -  Nogyi^a 
a  -  M^oynM  +  MpOYNP  +  MpovnP  +  Moowa  +  MooYHa 

where : 

Modyn’  M^dyni  Nodyn  Sind  are  defined  in  equations  44  to  47,  respectively,  and 
S^moyn  ~  Nrsin^oo  *  Sfpcos^oo  *  (Np  *  Sfr)sinaooosao  *  (Nqsinoo  *  Sf^sao)tanpo 

St^DYN  ■  ^Npsin^oo  -  Ucos^ooJ  t  ♦  (Lp  “  Nr)sinaoco8ao)secpo 
^aoYN  ”  (^qcosoo  *  Nqsinoo^secPo 

^  Note  this  is  a  simplified  equation  set  used  by  Kalviste  to  illustrate  computation  of  the  stability 
parameters.  The  complete  equation  set  is  given  in  Reference  (5). 


(56) 

(57) 

(58) 

(59) 

(60) 
(61) 
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M|4dyn  •  [  Mq  -  NrSin^ao  -  Sfpcos^ao  -  (Np  +  ^f)sinaocosao  (62) 

-  (^qcoscxo  *  Nqsinao)tanpo]sinPo  (MrSinao  Mpcosao)cospo 

MpoYN  “  MpSinao  -  Mrcoscxo  -  (Npsin^ao  -  Ljcos^oo  +  (Lp  -  Nr)sincxocosoo]tanpo  (63) 

MaovN  ‘  Mq  -  ^g'qcosao  NqsinaO)tanpo  (64) 

NpoYN  =  (  “  a!rSin^ao  +  Npcos^ao  +  (Nr  -  ^p)sinaocosoolcospo  (65) 

+  (Nqcosao  -  £fqSinao)sinPo 

N'pDYN  ”  ■  ifpsin^oo  “  Nrcos^ao  +  (Np  +  ^r)sinaocosao  (66) 

Nooyn  “  NqCOSao  -  SfqSinao  (67) 

Taking  the  Laplace  transform  of  equations  56  to  58  and  placing  them  in  matrix  form  yields, 


®  “  ifpDYN 

-  SfpDYNS  -  SfpoYN 

“  SfaovN®  “  ^^OOYN 

P' 

Njiqyn 

S^  +  NpDYNS  +  NpoYN 

Nodyn®  Nooyn 

P 

•0 

(68) 

-  MpQYN 

-  MpOYNS  -  MpoYN 

2 

S  ■■  Mqoyn^  ~  Modyn 

a 

» 

If  the  off-diagonal  terms  (coupling  terms)  of  the  above  matrix  equation  are  zero,  then  the 
diagonal  terms  represent  the  roll  subsidence,  Dutch-roll  and  pitch  short  period  stability  modes, 
respectively  (Reference  (5)).  Taking  advantage  of  the  fact  that  a  square  matrix  with  distinct 
eigenvalues  can  always  be  diagonalized  (^*  diag  [X.i,  Xa,  ...X.n])by  a  similarity  transforma¬ 
tion  (see  Reference  (24)),  Kalviste  rearranges  equation  (68)  as  shown  in  equation  (69)  with  the 
five  new  stability  parameters  (in  dimensional  form)  on  the  matrix  diagonal. 
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s-S£i 


MCOP 


0  S  +  Npcop 


0 

0 

2 

S  “  MocoP^  ~  ^OCOP 


McoP 

Pcop 

“cop 


(69) 


To  relate  the  new  stability  parameters  to  the  conventional  aircraft  modes  of  motion  Kalviste 
developed  a  method  of  integrating  the  characteristic  roots  (s>&-^j(o)from  the  uncoupled  equa¬ 
tions  to  the  coupled  equations.  To  account  for  coupling  between  the  conventional  modes  of  mo¬ 
tion,  multipliers  Ki  and  K2  are  placed  on  the  off-diagonal  coupling  terms  as  shown  in  equation 
(70). 


S  “  ii^MOYN 

“  SCpOYK*  ~  SfpDYN 

~  ^®OYN®  "  SfooYN 

K2(NHDyN) 

*  NpQYN*  NpQYN 

NaoYN®  NodyN 

p 

-0 

(70) 

K2(  Mjjqyn) 

Ki(~  MpoYN^  ~  ^Pdyn) 

2 

s  "  ModYN®  ”  ^ODYN 

a 

The  migration  of  the  characteristic  equation  roots  due  to  coupling  can  be  tracked  by  integrating 
the  characteristic  roots  with  respect  to  these  multipliers  from  zero  to  one  as  shown  in  Reference 
5.  Integrating  the  characteristic  roots  with  respect  to  Ki  accounts  for  coupling  between  the  lon¬ 
gitudinal  and  lateral/directional  modes  and  integrating  the  characteristic  roots  with  respect  to  K2 
accounts  for  coupling  between  the  roll  subsidence  and  the  Dutch-roll  modes,  in  this  manner  the 
migration  of  the  eigenvalues  can  be  tracked,  adding  insight  to  the  effects  of  lateral/directional  as 
well  as  longitudinal-lateral/directional  coupling. 
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The  coupled  stability  parameters  given  in  equation  (69)  were  derived  from  the  equations-of-mo- 
tion  that  included  the  linearized  aerodynamic  derivatives  (both  static  and  dynamic)  and  the 
kinematic  coupling  of  these  derivatives.  To  define  an  aircraft’s  aerodynamic  stability  for 
maneuvering  flight,  Kalviste  includes  the  inertial  effects  (see  equations  (71)  to  (73))  due  to 
steady  rotation  rates  in  the  equations-of-motion. 


Pi  -  I^QoR  +  RoQ] 


(71) 


Q; 


^Jjg=j[RoP.PoR] 


(72) 


Ri-  [PoQ  +  QoP] 


(73) 


To  accomplish  this,  the  appropriate  aerodynamic  derivatives  are  modified  to  account  for  the  iner¬ 
tial  coupling  effects  such  as; 


>  d 


TYPE  OF  CRITERION : 

The  five  stability  parameters  developed  by  Kalviste  relate  directly  to  the  stability  of  the  natural 
aircraft  stability  modes  (i.e.,  the  pitch  short  period,  roll  subsidence  and  Dutch-roll).  The 
parameters  can  thus  be  used  to  predict  open-loop  stability  of  an  aircraft  in  steady  maneuvering 
flight  (i.e.,  dynamic  aerodynamic  effects  and  kinematic  and  inertial  coupling  terms  are  included 
in  the  parameter  formulation).  Failure  to  satisfy  the  stability  criterion  of  one  of  the  five  stability 
parameters  has  not  been  correlated  with  aircraft  departure  susceptibility.  However,  regions  of 
the  a-p  plane  for  which  one  of  the  stability  criterion  is  not  satisfied  reveals  a  region  of  local 
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instability  which  results  in  degraded  flying  qualities  of  the  unaugmented  aircraft  and  possible 
departure  prone  characteristics. 

CRITERION  PLANE  : 

No  departure  susceptibility  criterion  has  been  correlated  with  the  Kalviste  Coupled  Stability 
Parameters.  Figure  46  shows  how  the  five  stability  parameters  can  be  plotted  on  an  a-p  plane 
to  determine  whether  any  regions  of  instability  exist. 

DESIGN  IMPACT : 

The  formation  of  new  dynamic  modes  due  to  the  effects  of  coupling  between  the  conventional 
aircraft  modes  can  be  identified  using  the  root  integration  technique  presented  by  Kalviste  in  Ref¬ 
erence  (5).  The  significance  of  this  method  is  that  it  applies  even  in  highly  coupled,  dynamic  con¬ 
ditions. 

CB1T£B1QN.SHQRTCQIVIINGS ; 

Kalviste's  formulation  of  the  equations-of-motion  is  such  that  the  effects  of  trajectory  stability 
(defined  in  terms  of  aircraft  translational  motion  due  to  aerodynamic  and  propulsive  forces  and 
the  force  of  gravity)  are  not  included.  This  excludes  the  possible  coupling  effects  between  the 
"fast"  modes  and  the  slow  trajectory  modes  consisting  of  the  conventional  phugoid,  spiral  and  al¬ 
titude  stability  modes.  Kalviste  admittedly  notes  that  inclusion  of  the  trajectory  effects  can  either 
stabilize  or  destabilize  the  aerodynamic  stability  characteristics  predicted  by  the  developed 
stability  parameters  of  equation  (69).  Kalviste  has  already  incorporated  the  trajectory  effects  in 
the  parameter  development  by  including  the  p,  y,  velocity  and  altitude  equations  in  the  formula¬ 
tion  of  the  equations-of-motion.  This  results  in  a  ninth  order  system  that  fully  encompasses  all 
coupling  effects  of  the  trajectory  into  the  already  developed  aerodynamic  stability  parameters. 

Levels  of  instability/stability  of  each  of  the  parameters  still  needs  to  be  further  investigated  to 
determine  whether  they  can  be  correlated  with  aircraft  departure  susceptibility. 
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Figure  46  Kalviste  Stability  Parameters  Plotted  Versus  Angle-of-Attack  and  Sideslip  Angle  (Reference  (5)) 
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CHODY:  ROUTH  CRITERION  PARAMETERS 
DATE/RESEARCHER: 

1989;  Chody  (Eidetics) 

CRJJERIQN: 

Extension  of  the  Cn^Qy^  and  LOOP  stability  concepts  to  address  asymmetric  (po  *  0;  but  (t>o  *  0) 
steady  maneuvering  flight.  Figure  47  illustrates  the  suggested  open-loop  stability  criteria  planes 
and  Figure  48  illustrates  the  suggested  dosed-toop  stability  criteria  that  refines  the  classical 
definition  of  the  LCDP  parameter. 

BASIS  OF  CRITERION: 

Motivated  by  inconsistent  correlation  of  recent  departure  experience  with  the  Cn^Qy^  stability 
parameter,  Chody  chose  to  linearize  the  lateral/directional  equations-of-motion  and  apply 
Routh's  Stability  Criteria  to  the  resulting  quartic  polynominal  charactaristic  equation  with  as  little 
simplification  as  possible  (see  equation  (21)).  As  was  shown  in  the  Cnp^yM  stability  criterion  sec¬ 
tion  the  expression  for  Cnp^y^  (see  equation  (22))  now  contains  two  additional  terms,  the  second 

of  which  (i.e.,  IC'jpC'n,  -  C'(,C'npl )  was  shown  in  Figure  22  to  be  significant  at  least  for  the 

case  of  the  preliminary  configuration  of  the  X-31 A  aircraft.  Note  that  equation  (22)  is  a  refine¬ 
ment  of  the  original  expression  derived  from  Cnpoy^  include  higher  order  terms  that  with 
today  and  future  aircraft  configuration  concepts  appear  to  be  important.  The  modified  expres¬ 
sion  for  CnpQy^  is  still  a  "Routh  Stability  Test”  on  the  s^  •  coeffident  (or  ”C”-coefficient;  see  equa¬ 
tion  (22))  of  the  uncoupled  lateral/diredional  diaracteristic  equation. 

In  addition  to  applying  Routh's  Stability  Criteria  to  the  ”C”  -  coefficient  to  determine  minimum 
levels  of  CnpQy^,  Chody  proposes  that  Routh’s  criteria  also  be  applied  to  the  other  coeffidents  of 
(he  lateral/diredional  charaderistic  equation  as  shown  in  Figure  47a  through  47c.  (Note,  Chody 
omits  addressing  the  Routh  discriminant,  0(BC-A0)-B^E;  Negative  values  of  the  discriminant 
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o- 

(a)  LateraVOiractional  Damping  for  Stability  Based  on  the  Characteristic  Equation  Coefficient 


(b)  Minimum  Value  of  Cojovn  Stabiiity  Based  on  the  s*  ("C^  Characteristic  Equation  Coefficient 


Figure  47  Open-Loop  Lateral/Directional  Stability  Boundaries  of  Chody  Based  on  the  Routh  Stability 
Criteria  (Reference  (33)) 
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ixx  Cnp 

Tu 


•  BOUNDARY  SHOWN  FOR  Ci,  ( I*  V  ♦  gooepo*!"*®)  *  ® 
•  •  BOUNDARY  SHOWN  FOR  Ci,(L',V-LVtm  So)  <  0 


(c)  Uteral-DirecJiofMJStabilityBoundariBsBasedon  s^fD" and -E^aiaracteristic Equation Ck)e(fid^^^ 


Figure  47  (Concluded).  Open-Loop  Laleral/Direcfional  Stability  Boundaries  of  Chody  Based  on  the 
Roulh  Stability  Criteria  (Reference  (33)) 
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Figure  48  Closed-Loop  Lateral/Directional  Stability  Boundaries  of  Chody  (Extension  of  LCDP  to  Include 
the  Aerodynamic  Dynamic  Deriivatives)  (Reference  (36)) 
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indicate  that  the  real  part  of  one  complex  pair  of  roots  (either  the  conventional  "Dutch-roll  pair  or 
the  lateral-phugoid")  is  in  the  unstable  RHP  which  implies  one  oscillatory  divergence  will  occur). 
Correlation  of  these  coefficient  stability  bounderies  with  aircraft  departure  susceptibility  has  yet 
to  be  done. 

The  same  analysis  approach  used  to  include  the  aerodynamic  dynamic  derivatives  in  the  expres¬ 
sion  for  CnpQy^  was  also  used  to  modify  the  LCDP  closed-loop  stability  parameter.  The  resulting 
expression  is  provided  in  equation  (30c)  and  is  plotted  versus  angle-of-attack  in  Figure  48  with 
some  simplification. 

TYRE  PP  CBITEBIQM: 

Refinements  of  the  open  and  closed-loop  lateral/directional  departure  susceptibility  parameters 
CnpQYN  LCOP  are  proposed.  Additional  open-loop  stability  bounderies  based  on  the  "B",  ”D" 
and  "E"  coefficients  of  the  iateral/directional  characteristic  equation  are  also  suggested. 

CRITERIQK£LQIS: 

See  Figures  47  and  48. 

DESIGN  IMPACT: 

As  suggested  in  Reference  (33),  "One  interpretation  of  the  modified  criteria  is  that  the  currently 
defined  minimum  value  requirements  for  CnpQYN  L.CDP  should  not  be  constant  for  all  con¬ 
figurations.  Instead  minimum  acceptable  values  should  be  based  on  the  configuration  specific 
dynamic  derivatives." 

The  proposed  criteria  is  an  important  refinement  because  it  is  consistent  with  modern  fighter/at¬ 
tack  aircraft  configuration  aerodynamics  which  are  dominated  by  the  effects  of  the  forebody  vor¬ 
tices  at  high  angles-of-attack  which  tends  to  result  in  adequate  static  stability  characteristics  but 
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poor  dynamic  stability  characteristics  (As  an  example  see  the  yaw  and  roll  damping  charac¬ 
teristics  for  the  X-31 A  given  in  Figure  4)  (Reference  (33)). 

Table  V  from  Reference  (33)  compares  the  major  differences  in  design  features  and  charac¬ 
teristics  of  the  fighter/attack  aircraft  of  the  1950’s  and  1960's  to  the  current  design  periods  begin¬ 
ning  with  the  early  1970’s  (i.e.,  F-14,  F-15,  F-16  generation). 

Review  of  the  substantiating  data  used  by  Weissman  (to  correlate  regions  of  versus 
LCDP  with  aircraft  departure  susceptibility)  (Reference  (33))  found  that  the  dynamic  data 
(Cnp.  C(p.  Cn,  and  C(,)  influenced  the  departure  susceptibility  characteristics  of  the  aircraft  in  cor¬ 
relation  data  base,  but  were  not  the  dominant  factors. 

Table  V  Comparison  of  the  Major  Differences  in  Design  Features  and  Characteristics  of  Fighter/Attack 
Aircraft  of  the  1950’s  and  1960's  versus  the  Current  Design  Period.  (Reference  (33)) 


50’s 'GO’S 

70'8  •  OO’s 

“Cw 

17^  23® 

30®.  40® 

Relatively  Low 

High 

Iz/lx 

Med  •  High 

High 

Forebody  Fineness  (I/d) 

2-2.5 

4-6 

Asymmetric  Yawing  Moment  (Cng) 

0-0.01 

0.05-0.12 

Directional  Weathercock 

Stability  (Cn#) 

Dominated  by  aft 
located  vertical  tail 

Dominated  by 
forebody  (nose/lex*/ 
canard)  vortioes 

Effect  of  Increasing  Cn^  on 
Directional  Damping  (Cn,) 

irtcreased  damping  due 
to  forces  acting  aft  of 
aircraft  CQ 

Decreased  damping 
due  to  forces  acting 
ahead  of  CG 

Impact  of  Rotary  Cross 

Derivatives  (C^.  Cnp) 

Low  due  to  relatively 

High  due  to  increased 

‘Leading  edge 
extension 
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As  experienced  during  early  X-31 A  configuration  design,  work  done  by  NASA  has  shown  that 
when  the  high  angle-of-attack  aerodynamics  of  a  configuration  are  dominated  by  the  forebody, 
there  is  a  fundamental  interchange  between  ttie  static  and  dynamic  stability  characteristics  (Ref¬ 
erence  (33)).  Specifically  if  the  forebody  is  altered  to  improve  the  static  directional  stability  (Cnp) 
at  high  angles-of-attack,  yaw  damping  (Cn,)  will  invariably  decrease.  Conversely,  aKerations  to 
the  forebody  to  improve  dynamic  stabiiity  characteristics  will  invariably  degrade  the  static  direc¬ 
tional  stability  characteristics. 

CRITERION  ADVANTAGES: 

The  process  of  determining  the  stability  characteristics  of  an  aircraft  without  computing  the 
roots  of  the  characteristic  equation  to  gain  insight  into  the  dominant  vehicle  characteristics  is 
a  formidible  analyticai  task  for  the  complete  linear  unaugmented  aircraft  equations-of-mo- 
tion.  To  make  this  a  realistic  task,  as  is  often  done,  Chody  utilizes  Kalviste’s  equation 
development  given  in  Reference  (51)  which  decouples  the  complete  linear  equations-of-mo- 
tion  and  considers  the  stability  characteristics  of  just  the  lateral/directional  set.  However,  un¬ 
like  the  original  CnpovN  LCDP  stability  parameters  derived  by  Moul  &  Paulson  (Reference 
(35)),  Chody  does  not  exclude  the  aerodynamic  dynamic  derivative  terms  (Cn„Cn„  etc.)  as 
second  order  effects  in  his  derived  expressions.  As  pointed  out,  the  influence  of  the 
dynamic  derivatives  can  be  significant  for  the  configurations  where  forebody  dynamics 
dominate  the  flow.  Furthermore,  by  application  of  Routh  Stability  Criteria  to  the  other  coeffi¬ 
cients^  of  the  lateral/directional  quartic  characteristic  equation  (A,  B,  C,  D,  E  and  D(BC-AD)- 
B^E  must  all  be  positive  for  absolute  system  stability),  Chody  has  found  that  dynamic 
instabilities  other  than  those  associated  with  the  "C"  -  coefficient  are  also  influential  to  an 
aircraft  configuration  being  departure  prone. 


^Note,  Chody  does  not  address  the  D(BC-AD)-B^E  term  of  Routh’s  Stability  Criteria 
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CRITERION  SHORTCOMINGS; 

Used  as  aircraft  departure  susceptibility  criteiia  the  lateral/directional  stability  criteria  proposed 
by  Chody  has  two  shortcomings  to  consider.  First,  as  Chody  suggests,  the  criteria  (with  the  ex¬ 
ception  of  the  extension  of  the  expressions  for  the  Cdm^vm  LCDP  parameters)  still  need  to  be 
further  analyzed  in  terms  of  correlating  the  proposed  parameters  with  aircraft  departure  suscep¬ 
tibility  characteristics. 

Secondly,  because  the  proposed  stability  parameters  are  derived  from  the  uncoupled 
lateral/directional  equation-of-motion  set,  instabilities  due  to  coupling  between  the  longitudinal 
and  lateral/directional  modes  of  motion  which  is  prevalent  at  high  angles-of-attack  can  not  be  ad¬ 
dressed.  If  significant  aerodynamic  cross-coupling  does  exist  (i.e.  Sfa.  Mp,  etc.)  the  criteria 
proposed  by  Johnston  of  STI  (^/TeJ  and  Kaiviste  (Npcop>^^oicof»  ^Ic.)  may  be  more  appropriate. 
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CONCLUSIONS 

The  goal  of  this  research  was  to  provide  guidelines  to  aid  in  the  analytical  determination  of  the 
departure  susceptibility  of  future  fighter  aircraft  designs.  Towards  this  end  the  major  departure 
susceptibility  criteria  developments  dating  back  to  the  concept  of  the  stability  derivative  were  dis¬ 
cussed  at  length  in  terms  of  their  derivation,  their  design  impact,  their  strong  points  and  their 
shortcomings.  As  expected,  no  one  analytical  criterion  can  be  applied  to  an  aircraft  design  and 
provide  the  "full  picture"  as  to  the  susceptibility  of  an  aircraft  to  depart  from  controlled  flight. 
However,  the  survey  of  criteria  revealed  that  although  in  some  cases  there  is  overlap  in  the  exist¬ 
ing  criteria,  on  the  whole  the  criteria  can  be  used  in  an  intergrated  fashion  (from  simple  to  more 
complex)  to  follow  the  design  evolution  and  the  associated  availability  and  refinements  of  the 
aerodynamic  da^a  base. 

Without  exception  each  of  the  departure  susceptibility  criteria  researched  were  based  on  linear 
or  quasi-linear  frozen  point  analysis.  Although  application  of  linear  analysis  does  still  appear  to 
apply  to  the  high  angle-of-attack  flight  regime,  caution  must  be  exercised  in  defining  limits  of 
validity  of  the  derived  linear  models.  Additionally,  where  aerodynamic  cross-coupling  exists 
(i.e.,  Cmp.  Cn,.  etc.)  and/or  inertial  coupling  is  an  important  stability  consideration,  the 
lateral/directional-longitudinal  equations-of-motion  can  not  be  decoupled,  and  the  linearization 
must  be  based  on  the  six  degree-of-freedom  equations-of-motion  and  trim  points  that  include 
asymmetric  flight  (Po  0, 4»o  0)  and  nonzero  angular  rates  (Po,  Qo,  Ro  *  0).  Furthermore,  the 
linear  analysis  based  results  must  address  both  open-  (characteristic  equation)  and  closed-loop 
(transfer  function  numerator  dynamics)  parameters  and  be  verified  using  nonlinear  simulations 
and  available  experimental  techniques.  Table  VI  summarizes  the  suggested  use  of  the  existing 
departure  susceptibility  criteria/stability  parameters  as  they  apply  in  the  course  of  the  design 
evolution. 
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Table  VI  Recommended  Methodology  For  EvaluatioiVDetermination  Of  Aircraft 
Departure  Susceptibility 


DATA/ANALYSIS 

REQUIREMENTS 

APPROPRIATE/RECOMMENDED  CRITERIA  AND  OR 
METHODOLOGY 

Inertial  Data 
(!»>  I2z>  Ixz) 

Determino  signHicanoe  of  aerodynamic  cross-coupling 
0.e.,  see  vector  polygon  method  of  Referenoe  (37 )). 

Lateral-Directional 
aerodynamic  static  data 

(C^,  Cftj,  Cnjj  . .  Cntf  CnJ 

a.  If  aerodynamic  cross-coupling  is  significant  apply 
KaMste's  Coupled  Stability  Criteria. 

b.  H  aerodynamic  cross-coupling  is  not  significant 
apply  STI/Weissman’s  LCDP  versus  Cn«QVN 

Criterion,  Pelikian’s  Cn^^  Criterion,  and  STI's  VT«,. 

Aerodynamic  dynamic  data 
(C^.  Cn,.  etc.) 

Apply  Chod/s  Stability  Criteria  to  extend  the 
concept  of  'Cnftfot*  to  include  the  effects  of  the 
dynamic  derivatives. 

Evaluate  the  effects  of 
asymmetric  flight  (ffo  *  0) 

Applicable  stability  parameters  include  STI’s  Vr«, 
and  KaMste’s  Coupled  Stability  Parameters. 

Assess/determine  the  effects 
of  maneuvering 

Kalvlste’s  (1989)  Dynamic  Stability  Parameters 
address  the  effects  of  steady  maneuvering  flight. 

Additionally  pilot-in-the-loop  simulation  is  highly 
recommended. 

Assess/determine  the 
effects  of  the  flight  control 
system  augmentation 

The  closed-loop  criteria,  VTt,  and  VT«,  .apply  equally 
well  to  fully  augmented  aircraft. 

Augmentation  effects  can  be  considered  in  Chody 
and  KaMste  fSS)  criteria  using  augmented  deriva¬ 
tives  (i.e.,  Lp^  •  L^  ♦  •  L^,  etc.).  Additionally 

the  LCDP  parameter  can  be  modified  to  include  the 
effects  of  an  aileron-rudder  interconnect. 

The  effects  of  the  FCS  ( especially  digital,  full- 
authority.  fly-by-wire  FCS)  are  best  assessed  using 
pilot-in-^e-loop  simulation- 
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RECOMMENDATIONS 


To  thoroughly  analyze  the  departure  tendencies  of  an  aircraft  configuration  (note,  even  small 
configuration  changes  can  alter  the  departure  characteristics  of  an  aircraft)  through  the  evolution 
of  the  design  (in  terms  of  data  availability  and  complexity)  each  of  the  following  factors  should 
be  considered: 


1 .  Analysis  of  both  open-  and  closed-loop  departure  susceptability  should  be 
addressed. 

RATIONALE:  Aircraft  departure  from  controlled  flight  can  occur  without  the 
pilot-in-the-loop  as  well  as  be  caused  by  pilot-in-the-loop  control.  If  closed-loop 
criteria  are  not  addressed,  aircraft  that  are  "assumed"  departure  resistant  can  in 
fact  be  just  the  opposite  when  the  effects  of  pilot  control  are  considered. 

2.  Careful  measurement  and  thorough  analysis  of  static  wind  tunnel  data 
should  be  conducted. 

a.  Consider  the  effects  of  Reynold's  number  during  high  angle-of-attack 
wind  tunnel  testing. 

b.  A  sufficient  number  of  data  points  must  be  taken  with  angle-of-attack  and 
sideslip  angle  to  reveal  nonlinearities  and  understand  the  limits  of  linear 
models  used  for  analysis  purposes. 

c.  Static  aerodynamic  cross-coupling  (i.e.,  Cmp,  Cio.Cna  etc.)  effects  should 
be  investigated  and  zero  sideslip  asymmetric  phenomena  (i.e.,  C^,  Cog) 
should  not  be  overlooked  or  ignored  as  anomallies. 

d.  Functionality  of  the  aerodynamic  data  should  be  carefully  analyzed,  e.g., 
is  control  data  a  function  of  sideslip  angle? 

RATIONALE:  Formulation  of  an  aerodynamic  math  model  that  accurately  rep¬ 
resents  the  aerodynamic  forces  and  moments  is  a  prerequisite  for  any  of  the 
proposed  criteria  to  successfully  correlate  with  aircraft  departure  susceptibility. 

3.  Include  the  influence  of  the  aerodynamic  dynamic  derivatives  in  the 
proposed  criteria. 

RATIONALE:  Current  and  future  fighter/attack  aircraft  designs  dominated  by 
forebody  vortices  at  high  angles-of-attack  can  lead  to  stable  statics 


« 


125 


NADC-90048-60 


(CnpQYN>())  unstable  dynamics  (Cn,  >0)  in  addition  to  the  dynamic  cross-cou¬ 
pling  derivatives  being  more  significant.  When  the  Cnp^yM  criteria  was 
developed,  aircraft  of  the  time  were  characterized  by  the  derivatives 
Yp,  Nr  andLp  being  negative  (stable).  The  opposKe  is  true  for  high  angles-of-at- 
tack  for  the  X-31 A  and  other  current  fighter/attack  aircraft  designs. 


4.  Conduct  a  sensitivity  analysis  of  the  high  angle-of-attack  aerodynamic 
characteristics  (i.e.,  Mp,  sfo.  Na,  Np,  Lp,  Ns,,  etc). 


RATIONALE:  The  reason  for  conducting  an  aerodynamic  sensitivity  analysis 
is  two  fold.  First  off,  the  dominant  aerodynamic  influences  on  the  aircraft’s 
departure  susceptibility  characteristics  can  be  determined  to  aid  in  addressing 
design  configuration  changes  and/or  establishing  a  flight  control  law  design  for 
departure  resistance.  Secondly,  the  sensitivity  analysis  provides  the 
aerodynamicist/f  light  controls  engineer  with  some  knowledge  as  to  the  robust¬ 
ness  properties  of  the  nominal  aerodynamic  math  model.  As  discussed  earlier, 
one  of  the  biggest  challenges  of  estimating  aircraft  dynamics  at  high  angle-of-at¬ 
tack  flight  conditions  is  the  development  of  an  "accurate"  aerodynamic  math 
model  that  can  be  used  with  a  high  degree  of  confidence.  The  sensitivity 
analysis  addresses  this  issue  and  highlights  where  model  uncertainty  is  critical. 


5.  Evaluate  the  effects  of  longitudinal-lateral/directional  dynamic  coupling  on 
aircraft  stability  which  are; 

a.  introduced  through  asymmetric  flight  conditions  (Po  ’‘0;  ^  *  0). 

b.  introduced  through  the  effects  of  maneuvering  (Po,Qo,  Ro  **  0;  i.e.,  inertial 
coupling  effects). 


RATIONALE:  Departure  susceptibility  criterion  parameters  such  as 
CnpoYN'Cnp*pp.  LCDP,  etc.,  based  on  decoupling  the  lateral/directional  equations- 
of-motion  from  the  longitudinal  set  are  based  on  linearizing  the  six  degree-of- 
freedom  equations-of-motion  about  a  trim  flight  condition  that  assumes  zero  roll 
and  yaw  angular  rates  as  well  as  zero  bank  angle  (typically  pitch  rate,  Qo  is 
also  assumed  to  be  zero). 


Typical  of  most  air  combat  maneuvers,  high  angular  rates  can  increase  an 
aircraft's  departure  susceptibility  from  two  perspectives.  Firstly,  the  pilot 
workload  is  increased  significantly  during  maneuvering  flight.  Secondly,  non¬ 
zero  pitch  rate  destabilizes  the  normal  modes  (primarily  dutch-roll)  and  reduces 
the  available  control  power,  while  nonzero  roll  rate  causes  lateral/directional-lon¬ 
gitudinal  coupling  which  produces  mode  shapes  unfamiliar  to  the  pilot. 
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6.  Evaluate  the  effects  of  any  flight  control  system  augmentation  using  the 
applicable  analytical  departure  susceptibility  criteria. 


RATIONALE:  In  light  of  tfte  fact  that  future  advanced  fighter/attack  airaaft  will 
invariably  contain  full  authority,  digital  fly-by*wire/light  augmentation  systems 
that  "guarantee"  stability  (provided  sufficient  control  power  is  available),  applica¬ 
tion  of  the  "static  augmented"  departure  susceptibility  criteria  (i.e.  Pelikan,  Kal- 
viste)  serves  to  optimize  the  unaugmented  aircraft  dynamics  including  the 
effects  of  the  controls  to  minimize  control  power  requirements  from  the  outset  of 
the  design  evolution. 


7.  Investigate  the  effect  of  nonlinearities  via  man-in-the-loop  simulation  to  in¬ 
clude: 

a.  the  effect  of  the  pilot  (i.e.,  determine  closed-loop  aircraft  pilot  departure 
susceptibility). 

b.  the  effect  of  nonlinearities  such  as  position  and  rate  limiters,  hysteresis, 
nonlinear  dynamic  maneuvering  effects,  atmospheric  disturbances. 

c.  evaluation  of  failure  states  and  degraded  flight  control  system  modes. 

d.  verifying/validating  the  linear  departure  susceptibility  crKeria  findings. 


RATIONALE:  Determining  local  regions  (in  terms  of  M,  h,  ao,  po)  of  stability/in¬ 
stability  associated  with  large  angle,  nonlinear  maneuvering  flight  can  become 
very  cumbersome  using  quasi-linear  methods  (the  problem  is  further  com¬ 
pounded  if  the  effects  of  the  flight  control  system  are  included).  In  addition, 
their  determination  might  not  be  indicative  of  aircraft  departure  susceptibility 
since  the  aircraft  trajectory  might  diverge  into  a  "stable  region"  where  the 
aircraft  becomes  controllable  again. 

Utilizing  a  high  fidelity  pilot-in-the-loop  simulation,  closed-loop  stability  of  the 
aircraft/pilot  system  can  be  determined  from  a  nonlinear  perspective.  "Flying"  a 
test  matrix  of  critical  departure  prone  maneuvers,  the  normal  aircraft  flight  trajec¬ 
tory  can  be  described  as  locally  stable  if  after  any  disturbance  (i.e.,  due  to  tur¬ 
bulence,  gusts,  system  failure,  etc.)  from  the  normal  trajectory  it  converges 
back  to  the  nominal  flight  path  and/or  stabilizes  at  a  new  trim  condition.  If  the 
aircrafVpilot  system  is  locally  unstable  at  some  point  along  a  desired  flight  path 
trajectory  then  the  aircraft  will  not  return  to  the  nominal  flight  trajectory.  If  the  in¬ 
stability  is  not  "controllable"  (i.e.,  the  aircraft  does  not  follow  conimanded  inputs) 
at  the  very  least  the  flying  qualities  would  be  degraded.  In  the  worst  case 
scenario  the  aircraft  would  diverge  from  the  desired  flight  trajectory  (and  not 
"stabilize"  in  a  different  trim  condition),  with  post-stall  gyrations  to  follow  and 
spin  susceptibility  very  high. 
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8.  Support  analytical  and  pilot  simulation  results  with  available  experimental 
methods. 


RATIONALE:  Because  of  the  limKations  of  winJ  tunnel  data  acquisition/inter¬ 
pretation  (Reynolds  number  effects,  amplitude/osdilation  effects,  etc.)  and  the 
inherent  limits  of  any  mathematical  model  { be  H  pilot-in-the-loop  or  otherwise), 
it  is  recommended  that  aircraft  departure  susceptibility  analysis  make  full  use  of 
experimental  techniques  such  as  water  tunnel  tests  (for  flow  visualization), 
tethered  model  tests  and  drop  model  testing  to  refine  and  validate  the  high 
angle-of-attack  and  departure  susceptibility  characteristics  of  the  subject  aircraft. 
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DEPARTURE  :  The  event  indicating  loss  of  control  which  may  develop  into  a  post-departure  gyration, 

spin  or  deep  stall  condition.  The  departure  may  be  characterized  by  divergent,  large 
amplitude,  uncommanded  aircraft  motions,  such  as  nose  slice  or  pitch-up.  An  AOA 
excursion  is  not  considered  a  departure  (Reference  (55)). 


DEPARTURE  SUSCEPTIBILITY  (MIL-S-83691 A  Definitions) ; 

Extremely  Susceptible  to  Departure:  Departure  from  controlled  flight  will  generally 
occur  with  the  normal  application  of  pitch  control  alone  or  with  small  roll,  and  yaw 
control  inputs. 

Susceptible  to  Departure:  Departure  from  controlled  flight  will  generally  occur  with 
the  application  or  brief  misapplication  of  pitch,  roll,  and  yaw  controls  that  may  be  an¬ 
ticipated  in  operational  use. 

Resistant  to  Departure:  Departure  from  controlled  flight  will  only  occur  with  a  large 
and  reasonably  sustained  misapplicaton  of  pitch,  roll,  and  yaw  controls. 

Extremely  Resistant  to  Departure:  Departure  from  controlled  flight  can  only  occur 
after  an  abrupt  and  inordinately  sustained  application  of  gross,  abnormal,  pro-depar¬ 
ture  controls. 


POST-STALL  GYRATION:  An  uncontrolled  oscillation  about  any  or  all  axis  following  a  departure. 


SEIN: 


A  sustained  rotation  in  yaw  at  an  angle-of-attack  greater  than  the  stall  angle-of-attack. 
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